Cascade atom in high-Q cavity: 
The spectrum for non-Markovian decay 



B.J. Dalton 1 and B.M. Garraway 2 

1 ARC Centre for Quantum-Atom Optics and Centre for Atom Optics and Ultrafast Spectroscopy, 
Swinburne University of Technology, Melbourne, Victoria 3122, Australia 
2 Department of Physics and Astronomy, University of Sussex, Falmer, Brighton, BN1 9QH, United Kingdom 

(Dated: February 1, 2008) 

The spontaneous emission spectrum for a three level cascade configuration atom in a single mode 
high-Q cavity coupled to a zero temperature reservoir of continuum external modes is determined 
from the atom-cavity mode master equation using the quantum regression theorem. Initially the 
atom is in its upper state and the cavity mode empty of photons. Following Glauber, the spectrum 
is defined via the response of a detector atom. Spectra are calculated for the detector located inside 
the cavity (case A), outside the cavity end mirror (Case B — end emission), or placed for emission 
out the side of the cavity (Case C). The spectra for case A and case B are found to be essentially the 
same. In all the cases the predicted lineshapes are free of instrumental effects and only due to cavity 
decay. Spectra are presented for intermediate and strong coupling regime situations (where both 
atomic transitions are resonant with the cavity frequency), for cases of non-zero cavity detuning, 
and for cases where the two atomic transition frequencies differ. The spectral features for Cases 
B(A) and C are qualitatively similar, with six spectral peaks for resonance cases and eight for 
detuned cases. These general features of the spectra can be understood via the dressed atom model. 
However, Case B and C spectra differ in detail, with the latter exhibiting a deep spectral hole at 
the cavity frequency due to quantum interference effects. 



I. INTRODUCTION 

Since the original paper of Purcell [l[ the modification 
of atomic radiative decay due to the presence of optical 
and microwave cavities, photonic band gap (PBG) mate- 
rials has become a familiar topic of study in quantum op- 
tics. The presence of these systems structures the modes 
of the electromagnetic (EM) field, resulting in situations 
where the mode frequency dependences of the mode den- 
sity or the atom-field coupling constants (or both) are 
changed from the free space situation. The quantum EM 
field acts as a structured reservoir whose behaviour de- 
pends on the product of the mode density with the square 
of the coupling constants — the reservoir structure func- 
tions. If the reservoir structure functions vary slowly 
with mode frequency, the corresponding reservoir corre- 
lation time may be short compared to the time scale over 
which the atomic state changes, and the atomic density 
operator satisfies a Markovian master equation. This 
results in the irreversible decay of excited atomic state 
populations, though with a decay rate different to that 
for radiative decay in free space. This regime is charac- 
terised by a weak coupling between the atom and the EM 
field. On the other hand, if the reservoir structure func- 
tions vary rapidly with mode frequency, the correspond- 
ing reservoir correlation time may be long compared to 
the time scale over which the atomic state changes, and 
the atomic density operator no longer satisfies a Marko- 
vian master equation. A reversible oscillatory decay of 
excited atomic state populations results, with a period 
related to the atom-field coupling constant. This regime 
is characterised by a strong coupling between the atom 
and the EM field. 

There are a large variety of effects associated with 



atoms in structured EM fields, including modifications to 
the spontaneous emission, fluorescence and probe absorp- 
tion spectra, and comprehensive surveys may be found in 
several reviews, for example @, d, El- Amongst the many 
effects is the modification of the spectrum of the radia- 
tion emitted when an excited atom decays into initially 
empty field modes. It has been known from the early 
days of quantum optics that for an initially excited two- 
level atom (TLA) in free space, this so-called spontaneous 
emission (SE) spectrum is a Lorentzian centred around 
the atomic transition frequency and with a width given 
by the excited atom decay rate § . In the weak coupling 
regime the spectrum will remain Lorentzian, but with 
a modified width reflecting the changed atomic decay 
rate. The situation changes qualitatively for the strong 
coupling regime, with the spectrum no longer remaining 
Lorentzian. The essential physics for a TLA coupled to 
a single mode cavity in the strong coupling regime can 
be described via the Jaynes-Cummings model [6], where 
the energy eigenstates allowing for the atom-cavity mode 
coupling are the so-called dressed atom states Q. 

The first calculation for the SE spectrum for a TLA 
in an ideal lossless cavity was carried out by Sanchez- 
Mondragon et al. Using a Heisenberg equation of 
motion approach and determining the spectrum from the 
two-time correlation function for the atomic dipole oper- 
ator, they found a double peaked spectrum with a sepa- 
ration between peaks given by the so-called one-photon 
Rabi frequency, that is the Rabi frequency associated 
with one photon in the cavity mode. Line shapes were 
associated with the spectrometer bandwidth Q , since no 
cavity bandwidth was included. Another early calcula- 
tion for a TLA in a high-Q cavity was that of Agar- 
wal who used the dressed atom model to predict 
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the SE and weak probe absorption spectra. This arti- 
cle determines the SE spectrum for different positions 
of the spectrometer, the spectrum associated with side- 
ways emission via direct SE being given in terms of the 
two-time dipole correlation function, whilst the spectrum 
associated with end emission via the cavity output mir- 
ror is given in terms of the two-time correlation func- 
tion for cavity creation and annihilation operators. Al- 
though qualitatively similar in exhibiting vacuum Rabi 
frequency splitting, the spectra are quantitatively differ- 
ent. The SE spectra shown had lineshapes determined 
for the situation where the spectrometer bandwidth was 
large compared to the cavity decay rate. Non-Lorentzian 
SE spectra were also found by Lewenstein et al. [ll[ for 
emission from a TLA into both high Q cavities and PBG 
systems, here an essential states approach was used and 
the spectrum defined in terms of the long time proba- 
bility for finding one photon in a cavity or background 
mode. U sing a master equation and quantum regression 
theorem [l2| approach, Carmichael et al. [l3| took into 
account both cavity damping and spontaneous emission 
damping to determine the side emission SE spectrum for 
a TLA in a high-Q cavity from the two-time dipole corre- 
lation function, based on a zero bandwidth spectrometer. 
Further discussion and references on the theory of the 
SE spectrum for a TLA in a high-Q cavity are given by 
Childs et al. [14j and Carmichael et al. [15| in the review 

No direct experimental measurements of the SE spec- 
trum for a single TLA in a high-Q cavity appear to have 
been carried out. However there are related experiments 
in the strong coupling regime, such as on the absorption 
spectrum for a weak probe field that demonstrates sin- 
gle atom vacuum Rabi splitting [ItJ , or single atom 
Rabi oscillations [HI], or the collapse and revival exper- 
iments with a single TLA that demonstrate cavity 
mode quantization [2(| Hl| . All these measurements can 
be interpreted in terms of the Jaynes-Cummings model 
Q and the associated dressed atom states 0j- The di- 
rect measurement of SE spectra would probably only be 
possible in the optical regime, since SE in the microwave 
regime would be too weak to detect. For the case of a sin- 
gle TLA in a high-Q cavity or a PBG system which is also 
coupled to a strong laser field, there is an extensive liter- 
ature dealing with the related spectra, in particular the 
spectra associated with resonance fluorescence or probe 
absorption. The work of Mollow [13] for the free space 
situation predicted a three-peaked (AC Stark effect) flu- 
orescence spectrum, and a dispersion-like probe absorp- 
tion spectrum. For a TLA in an ideal lossless cavity in a 
coherent state Agarwal et al. [23( predicted a large num- 
ber of spectral lines in the fluorescence spectra, reflecting 
the contributions from the numerous dressed atom states 
involved. These spectra can also be interpreted via the 
dressed atom model, see the recent reviews @, [3, 0] . 

In contrast to the case of the SE spectrum for a TLA 
in a high-Q cavity or a PBG system, the case of three 
level atoms (3LA) has received little attention. The SE 



spectrum for a 3LA in a lambda configuration in a PBG 
system has been treated by John et al. [24j via the essen- 
tial states approach for various detunings from the band 
gap edge, the spectrum (which shows doublet and win- 
dow effects) being defined via the photon emission prob- 
ability. Ashraf [25| determined the SE spectrum for a 
3LA in a lambda configuration in an ideal lossless cavity 
using the dressed states approach, the spectrum (which 
shows a doublet structure reflecting the one photon Rabi 
frequency) being given via the two-time dipole correla- 
tion function. Line shapes were due to the spectrometer 
bandwidth. For a 3LA in a cascade configuration in a 
PBG system, Bay et al. [26| used the essential states ap- 
proach to determine the SE spectrum for various detun- 
ings of one of the atomic transitions from the band gap 
edge, the latter (which is strongly non-Lorentzian) being 
defined via the photon emission probability. The treat- 
ment applies when the two atomic frequencies are quite 
distinct, enabling only one at a time to interact strongly 
with the PBG. Paspalakis et al. [27| treated the case of 
SE from a 3LA in a lambda configuration in a PBG sys- 
tem, using the same approach and spectrum definition 
as [24l | . Similar spectral features were found, including 
the window not previously noted by John et al. The 
case of a 3LA in a cascade configuration in an ideal loss- 
less cavity has been studied by Zhou et al. [28[ using a 
dressed atom approach and both end and side SE spectra 
were determined. Line shapes were due to the spectrom- 
eter bandwidth, as no cavity damping was included. The 
cavity mode was resonant with the average of the atomic 
transition frequencies. For the situation where the two 
atomic transition frequencies were equal, six peaks were 
found in both the side emission and end emission SE spec- 
tra, though in the latter case two of these were negligi- 
ble. When the two atomic transition frequencies differed, 
there were eight peaks in the side emission spectrum. 

In a recent paper [29| we have considered the non- 
Markovian decay of a three level cascade atom with both 
transitions coupled to a single structured reservoir of 
quantized field modes for both a high-Q cavity and in 
a PBG system. Based on the approach given in [3(| 
the dynamics of this system has been treated via the 
essential states approach, using Laplace transform meth- 
ods applied to the coupled amplitude equations. Non- 
Markovian behaviour for the population dynamics of the 
atomic system was found, such as oscillatory decay for 
the high-Q cavity case and population trapping for the 
photonic band-gap case. A Markovian master equation 
approach was also applied, in which the atomic system 
was augmented by a small number of discrete quasi- 
modes or pseudomodes, which in the quasimode treat- 
ment themselves undergo Markovian relaxation into a flat 
reservoir of continuum quasimodes. For the high-Q cav- 
ity case a single discrete quasimode was involved, for the 
PBG case two coupled discrete quasimodes were needed. 
The essential states and Markoff methods gave identi- 
cal results, showing that complicated non-Markovian be- 
haviour can be treated by enlarging the non-Markovian 
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system, thereby turning a non-Markovian problem into a 
Markovian one. 

In the present paper we now consider the SE spec- 
trum for the case of a cascade atom in a high-Q cav- 
ity. The spectrum is much richer than for the TLA case 
and quantum interference effects may now occur. Quan- 
tum coherence and interference phenomena are central 
in many applications of fundamental quantum optics re- 
sults, and a comprehensive review of such effects may 
be found in a special issue of this journal on quantum 
inteference 31]. Unlike previous work, we will consider 
the case of an ideal spectrometer with zero bandwidth 
in order to display the actual SE spectral line shapes 
without these being masked by instrumental effects. A 
master equation approach incorporating cavity decay will 
be used to evaluate the spectra via the quantum regres- 
sion theorem. Both end (Case A(B)) and side (Case C) 
SE spectra will be calculated, and we aim to exhibit the 
expected quantum interference effects that can occur in 
a cascade system. We will show that the spectrum in 
both cases are given in terms of the Laplace transforms 
of both the atom-cavity mode density matrix elements 
and the evolution operator matrix elements. 

Before developing the mathematical formalism it may 
be useful to consider the physical processes that can be 
involved in the decay of the cascade atom in a high-Q cav- 
ity and the registration of a photon arrival in a suitable 
detector — this will enable possible interference effects to 
be identified. In a cascade atom initially in its upper 
state 1 2) and coupled to a high Q cavity empty of pho- 
tons, two photons can be emitted before the atom makes 
a transition to its lowest state |0) via the intermediate 
state |1). If there is a detector atom weakly coupled to 
the atom-cavity mode system and prepared in its lower 
state \A), the detector atom could absorb one of these 
photons and make a transition to its upper state \B). 
Following the approach of Glauber [32] the overall tran- 



sition probability for the \A) — > \B) process considered 
as a function of the detector atom transition frequency lo 
may be used for an operational definition of the sponta- 
neous emission spectrum for this cascade atom in a high 
Q cavity — the cavity mode being coupled to empty ex- 
ternal modes via the cavity mirror. Since the detector 
atom is only weakly coupled we need only consider pro- 
cesses with a single \A) — > \B) transition. Also, steps 
where a photon is created in an external mode would 
not be reversible. On the other hand if the atom-cavity 
mode coupling is strong, atomic transitions accompanied 
by photon number changes in the cavity mode may be 
reversible. If we consider the case {Case A) where the 
detector atom is placed inside the cavity, then the pho- 
ton causing the detector atom transition must have come 
from the cavity mode. However photons in the cavity 
mode can also be transferred to an external mode via loss 
through the cavity mirror, and if one photon is used to 
cause the \A) — > \B) transition in the detector atom, the 
other will be transferred to an external mode. We con- 
sider states of the combined detector atom, cavity mode, 
cascade atom, external mode system of the product form 
\D; n; v; m), where D — A, B; n = 0, 1, 2, . . .; v = 0,1,2; 
m = 0, 1, 2, . . . specify the detector atom state, the cavity 
photon number, the atomic state and the photon number 
for a specific external mode respectively. Then the over- 
all process in which the cascade atom changes from upper 
state 1 2) to lowest state |0), the detector atom changes 
from lower state \A) to upper state \B) and one pho- 
ton appears in a specific external mode, whilst the cavity 
mode is initially and finally empty of photons is denoted 
\A; 0; 2; 0) — » \B; 0; 0; 1). However this overall process in- 
volving the overall emission of two photons (one being 
absorbed by the detector atom, the other appearing in 
an external mode) has a number of different quantum 
pathways: 



|v4;0;2;0) «-> \A; 1; 1; 0) ~~> \A; 0; 1; 1) «-► \A; 1; 0; 1) — > \B; 0; 0; 1) 

|v4;0;2;0) «-> \A; 1; 1; 0) -f \B; 0; 1; 0) <- \B; 1; 0; 0) * \B; 0; 0; 1) 

|v4;0;2;0) «-> \A; 1; 1; 0) «- \A; 2; 0; 0) ~» | A; 1; 0; 1) -+ \B; 0; 0; 1) 

|v4;0;2;0> «-> \A; 1; 1; 0) ~ \A; 2; 0; 0) \B; 1; 0; 0) ~» \B; 0; 0; 1) (1) 

Reversible atom-cavity mode transitions are designated <->, irreversible cavity-external mode transitions ~~> and weak 
detector atom transitions — >. In addition, each of the reversible steps in the above pathways may involve further 
pathways if the atom-cavity mode coupling is strong. For example, the first pathway 

\A; 0; 2; 0) «-> \A; 1; 1; 0) - \A; 0; 1; 1) ~ \A; 1; 0; 1) \B; 0; 0; 1) (2) 

could branch into two different sub-pathways 

\A; 0; 2; 0) <-> \A; 1; 1; 0) «-> \A; 2; 0; 0) <- \A; 1; 1; 0) - \A; 0; 1; 1) <- \A; 1; 0; 1) -> \B; 0; 0; 1) (3) 
| A; 0; 2; 0) «-> \A; 1; 1; 0) «-> \A; 0; 2; 0) ~ \A; 1; 1; 0) - \A; 0; 1; 1) ^ \A; 1; 0; 1) -> \B; 0; 0; 1) (4) 

I 



The possibilities for branching are endless. The tran- sition amplitudes for all these different pathways (and 
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sub-pathways) are added to form the overall transition 
amplitude, and hence the transition probability would be 
expected to exhibit quantum interference effects. Thus 
the spontaneous emission spectrum for a three level cas- 
cade atom should demonstrate interesting interference 
phenomena. 

Of course there will also be a detection processes in 
which the cascade atom changes from upper state |2) to 
intermediate state |1), the detector atom changes from 
lower state \A) to upper state \B) but no photon appears 

I 



if the atom-cavity mode coupling is strong . For this 
overall process there still may be quantum interference 
effects. Similar processes to these would apply if the 
only atomic states considered were |2) and |1), with the 
coupling to the lower state |0) being set to zero, and sub- 
pathways involving |0) excluded. In all these cases the 
amplitudes for different sub-pathways would combine to 
produce an overall transition amplitude, and hence inter- 
ference effects giving maxima and minima in the detec- 
tion probability (spectrum) could occur. Whether these 
interference effects can be be demonstrated directly via 
altering the system coupling constants, detunings etc. is 
a separate issue. 

A similar discussion can be presented for the case 
where the detector atom is placed just outside the cav- 
ity mirror (Case B) to detect end emission, or where the 
detector atom is placed outside the cavity to detect the 
weak side emission (Case C). In case B it is the exter- 
nal mode photon absorption that is associated with the 
detector atom transition, in case C it is the transition in 
the cascade atom itself, so the details of the discussion 
will differ for these two cases. 

In Section |TT] the general features of the quasimode 
description used to treat the cascade atom-cavity mode 
system are covered, along with a brief outline of the ap- 
proach for defining the spectrum and how the quantum 
regression theorem is used in the calculations. Section liTEl 
covers the master equation for the cascade atom-cavity 
mode system, the related evolution operator and their 
determination in terms of Laplace transforms. Expres- 
sions for the SE spectra for three different positions of 
the spectrometer atom — inside the cavity (Case A), out- 
side the cavity to detect end emission (Case B) and out- 
side the cavity to detect side emission (Case C) are given 
in Section IIIIl Numerical results for the SE spectrum 
for the cascade atom in a high-Q cavity are presented in 
Section IIVI and conclusions given in Section [V] Details 
are presented in the Appendices. 



in a specific external mode, whilst the cavity mode is 
initially and finally empty of photons. This is denoted 
|A;0;2;0) — > |£?;0;1;0) and involves only one photon 
emission, the emitted photon being absorbed by the de- 
tector atom. Here 

|v4;0;2;0) *-> |A;1;1;0) -> |B;0;1;0) (5) 

is the only quantum pathway, though again there are 
numerous sub-pathways such as 
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II. GENERAL THEORY 

A. Atom-quasimode system 

The general system of interest is a radiating atom with 
energy states \E a ) coupled to a discrete set of quasimodes 
i of the quantum EM field, which in turn are coupled 
to a continuum set of quasimodes A. The Hamiltonian 
H for the atom-quasimodes system is given as the sum 
of the atomic Hamiltonian Ha, the quasimodes Hamilto- 
nian Hq and an interaction between the atom and quasi- 
modes Haq H3| 

H = H A +H Q +H AQ , (7) 

where 

H A = J2^\E a )(E a \ (8) 

a. 

Hq = Hui ateti + ^2 ftVij a\a,j 

+ J dA Pc (A)hAb(Ayb(A) 

+ J2 JdAp c (A) (hWi(A)atb(A) + H.C.) 

■i 

(9) 

= E \E a )(Ef,\ + H.C.). (10) 

E a >E{3 i 

In these equations E a = hui a is the atomic energy, Vi 
and A are the frequencies of the discrete and continuum 
quasimodes, aj,ai and 6(A)t,6(A) are the standard cre- 
ation, annihilation operators for these quasimodes with 
commutation rules [ai,aj] = [6(A), 6(A') t ] = <5(A — 
A')/pc(A), hVij is the coupling energy between discrete 
quasimodes, hWi(A) describes the coupling between dis- 
crete and continuum quasimodes and h\i- a /3 specifies the 



L4;0;2;0) <-► \A; 1; 1; 0) <-► \A; 2; 0; 0) ~ \A; 1; 1; 0) \B; 0; 1; 0) 
|A;0;2;0) «-> L4; 1; 1; 0) <- \A; 0; 2; 0) <- \A; 1; 1; 0) -f \B; 0; 1; 0) 



coupling between the discrete quasimode and the atomic 
transition between states \E a ) , \Ep). All interactions are 
given via the rotating wave approximation. The contin- 
uum mode density is pc{A). 

In the case where Wi(A) and pc{A) are slowly vary- 
ing functions of the continuum quasimode frequency, the 
system consisting of the atom plus discrete quasimodes 
is described by a reduced density operator p satisfying a 
Markovian master equation (30| 

d i — 

g- t P = ~j[Hs,p\ +22^PcW i W*([a j ,p4] + [ajp, a}}), 

ij 

(11) 

where the first term involves the Hamiltonian for the 
atom-discrete quasimode system H$ and the second term 
describes relaxation due to the interaction between the 
discrete quasimodes and a reservoir consisting of the con- 
tinuum quasimodes. The Hamiltonian H$ is given by 

H s = Ha + H qd (12) 
Hqd = hvi alfii + hVii a%a>i- (13) 

B. Spectrum 

The spectrum S(u>) is defined operationally [H, [33| in 
terms of the long time transition probability to the upper 
state \B) for a two- level atom detector with transition 
frequency to that is initially in its lower state \A) and 
which is weakly coupled to the EM field, and where the 
final state of the atom-quasimode system is unobserved. 
The Hamiltonian H for the combined system of detector 
atom and the radiating atom-quasimode system is then 

H = H + H D + V SD (14) 

where the Hamiltonians for the detector atom and its 
coupling to the radiating atom-quasimode system are 

H D = HujSz (15) 
V SD = KV+S- + V-S+). (16) 

Here S+ = \B) (A\, S- = \A) (B\ and S z = §(|B) (B\ - 
\A) (A\) are the usual atomic spin operators. The op- 
erators V+, V- depend on where the detector atom is 
situated and hence what region of the quantum EM field 
it samples. In all cases V+, V- are proportional to the 
negative, positive frequency components £+ of the 
electric field operator at the detector atom. For the case 
where the detector atom is situated inside a high Q cav- 
ity or is situated outside the output mirror, the EM field 
can be conveniently described in terms of quasimodes 
[13, [34, 35]. In the case {Case A) when it is inside a high 
Q cavity — see figure Q] — the detector atom samples the 
discrete quasimodes, in another case (Case B) when it is 
situated outside a high Q cavity to detect end emission — 
see figure Q] — it samples the continuum quasimodes. For 
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Case A Case B Case C 



FIG. 1: The cascade atom in a high-Q cavity with the detec- 
tor atom in various locations. In Case A the detector atom 
is inside the cavity, in Case B it is outside the cavity output 
mirror and positioned to detect end emission and in Case C 
it is outside the cavity positioned to detect side emission. 

cases when the detector atom is situated outside a high Q 
cavity to detect side emission — see figure Q] {Case C) — 
the EM field is conveniently described in terms of true 
modes [H, [37} • In the three cases we have 

V- = ^/i|ai = (7 + ) f Case A (17) 

i 

V- = J dA Pc {A)p*{A)b{A) = {V+Y 

Case B (18) 
V- = ]T K P \E a ) (E P \ - (?+)t 

E a <Ep 

CaseC, (19) 

where pi, p{A) and R a p are weak coupling constants. 
For simplicity we will ignore retardation effects — these 
could be included by incorporating the position of the 
detector atom in p*, p*{A) and R a p and that of the 
radiating atom in the X*. a g, and by arranging the de- 
tector atom to begin responding when the emitted EM 
field first reaches it. In Case C the expression for V- 
is based on the long distance radiating atom contribu- 
tion to the EM field positive frequency component £+ 
being related to the electric dipole operator of the ra- 
diating atom (ignoring the free evolution term, which 
gives a zero contribution to the spectrum), and V— can 
be expressed in terms of the downward atomic transition 
operators \E a ) {Ep\ {E a < E fi ). 

If the initial state of the atom-quasimodes system is 
a pure state |5>,) and the detector is in state \A) then 
the probability amplitude for the transition to an atom- 
quasimode system state |$/) and the detector in state 
\B) is given by 

A FI = (F\ U{t) \I) , (20) 

with \I) = |$,) \A) and \F) = |$/) \B) and 
U(t) — exp{—iT-Lt/h) the evolution operator for the com- 
bined atom-quasimode-detector system. If the atom- 
quasimode system is initially in a mixed state with den- 
sity operator jo} = J2Pi (^1 then the transition 

i 
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probability to the upper detector atom state and to any 
final state of the atom-quasimode system is given by 

^A = ££pi|<iW)|J>| 2 , (21) 

/ * 

and for long times t this will be taken as the un- 
I 



(F\U(t)\I) 



and finally as t — > oo 

S{uj) = ^2 J J dtidt 2 expiuj(ti - t 2 ) 

xTr (piV+(t 2 ) £_(<!)), (23) 

where 

VL(t) = exp(+iHt/h)V- exp(-iHt)/K) 

V+(t) = exp(+iHt/h)V+ exp(-iHt)/K) 

= (V-(t)) f (24) 

are the system Heisenberg picture operators associated 
with V- and V+ . The trace is over the radiating atom- 
quasimodes system. The expression for the spectrum is 
real, as expected. This approach defines the so-called 
operational spectrum. 

A rather different approach to defining the spectrum 
is based on the essential states approximation. Here 
the state vector for the full atom and EM field modes 
system — the modes are still described via quasi-modes — 
is expanded in a basis set \E a ; m; m(A)) where rii and 
m(A) are the numbers of photons in the ith discrete and 
A continuum quasimodes. Thus 

!*(*)) = EE / dA /MA) KE a ;n t ;m(A);t) 

a rii m(A) 

x exp(— i{uj a + riiVi + m(A)A}i) 
x\E a ; ni ;m(A)), (25) 

with interaction picture amplitudes b(E a ;rii;m(A);t). 
The long time one photon emission spectrum would be 
defined in this approach via the long time probability of 
finding one photon in a continuum mode of frequency A 

S/(A)=^^|&(£ Q ;^;l(A);oo)| 2 , (26) 

in which the continuum quasimode frequency acts as 
the spectral variable. This approach may be referred 
to as defining the ideal spectrum, corresponding to quan- 
tities arising naturally from fundamental considerations 



normalised spectrum S(u>). 



A straightforward quantum treatment [3a ] which in- 
volves expressing the evolution operator as a sum of 
terms with increasing orders of Vsd then gives 
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of what can in principle be measured in quantum me- 
chanics, rather than to the analysis of the behaviour of a 
model spectrometer. An alternative definition of an ideal 
spectrum [39| could be based on true mode states — which 
also form a continuum — rather than continuum quasi- 
modes states. The interrelationship between the differ- 
ent approaches to defining spectra has been considered 
by Cresser [33| . Calculations of the ideal spectrum would 
involve solving the coupled amplitude equations. 

C. Quantum regression theorem 

Determination of the operational spectrum requires 
the evaluation of two-time correlation functions, and here 
the quantum regression theorem is used. The usual state- 
ment of the quantum regression theorem [l2|, E3, EH is 
that for a system operator Yi(t) in the Heisenberg picture 
whose average satisfies the linear equation (t 0) 

|(%)) = E G ^')(^(*))< ( 2? ) 

3 

with matrix elements Gij (t) , then we can assert that the 
two time correlation functions (r ^ 0) satisfy a similar 
linear equation involving the same matrix elements 

(28) 

j-^mit + r)) = J2G^(r)(Yim(t + r)). 

(29) 

A useful result giving the two time correlation func- 
tion (^Sfj a (t + T)Ss~f(t)^ for system transition operators 

in terms of matrix elements of the evolution operator and 
density matrix elements [4l| can be derived from this 
form of the quantum regression theorem by considering 
the case where Yi is a system transition operator Yi = 



= — exp(-iu)t) / dhexpHuh) ($ f \exp(-iH(t- tJ/hW-expt-iHCtd/h) $j) , (22) 
ifr Jo 
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Sp a = \0) (a\. In this case (9i(t)^ = (Sp*(t)) = p a p[t), 
the density matrix element which satisfies the master 
equation. In general we can write the solution to the 
master equation in terms of matrix elements U a p- jS (i) 
of the evolution operator 



Pa/3 

7<5 



(30) 



We then can show that for (t ^ 0, r 0) the two impor- 
tant results 



S/3ct(* + T)S« 7 (tn = ^ Ugf3;5n (t) Ayjj(t) (31) 

S 0a (t)S 5l (t + T)\ = ^2 U fS;^{T)p^{t). (32) 



These results give the two time correlation function 

(Sp a (t + r)Ssy(t)^ for system transition operators in 

terms of matrix elements of the evolution operator and 
density matrix elements. 



III. THREE LEVEL CASCADE ATOM IN 
SINGLE MODE CAVITY 

A. Master equation 

The master equation for a three level cascade atom in 
a single mode high Q cavity coupled to the reservoir of 
continuum quasimodes at zero temperature is given by 

^p = ~[H s ,p[ + lr([a,pa^} + [ap, a% (33) 

where the Hamiltonian for the three level atom plus cav- 
ity mode system in the rotating wave approximation is 



hLo c a'a + h(u)Q — e>)(722 — ft(^o + <5)<5oo 
+h\a) (52S7 +gi&i) + {92S2 +9i&f) a]- 



(34) 

Here a , a' are the mode annihilation, creation opera- 
tors for the single cavity quasimode with frequency u> c . 
The operators a+ = |2) (1|, aj = |1) (2|, d+ + = |1) (0|, 
&i_ = |0) (1| and (T22, an, ffoo are atomic transition and 
population operators respectively, involving the upper 
state 1 2), the intermediate state |1) and the lower state 
|0). The lower (1 <-> 0) and upper (2 <-> 1) atomic tran- 
sition frequencies are uj\, L02 respectively. The average 
atomic transition frequency is ujo and S is half the dif- 
ference between lower and upper transition frequencies. 
Thus ujq = (iui + W2V2, 5 = (uji — u>2)/2. The atomic 
term in the Hamiltonian is based on adding a constant 
^(5 : 22 + <7n + f?oo) to the original Hamiltonian so as to 
make the energy of the intermediate state |1) equal to 
zero. The atomic states are illustrated in figure [H along 



hco, 



12) 



hco, 



ID 



hco 



10) 



IB) 



IA> 



FIG. 2: Energy levels and energy differences in the three 
level cascade atom (left) . The offset 5 represents the difference 
between the energy of level 1 and the midpoint between levels 
and 2. The energy levels for the two state detector atom 
are also shown (right). 



with the states of the detector atom. The one photon 
Rabi frequencies are 52 and g\ for the two transitions 
and r is the cavity decay rate, where we have used the 
notation A21 — > 32, A10 — > gi- The coupling constants g\ 
and g 2 are proportional to the scalar product of the vector 
dipole matrix elements (2| |1), (1| |0) between the 
upper and intermediate states or the intermediate and 
lower states with the polarization unit vector for the cav- 
ity quasimode, and by a suitable choice of phase can be 
taken as real and positive. The reduced density operator 
is p. In terms of the flat continuum mode density pc and 
discrete-continuum quasimodes coupling constant W, we 
have T — 2tt pc \ W\ 2 . At zero temperature the contin- 
uum quasimodes are empty of photons. In the present 
case the atomic spontaneous emission rate 7 into side- 
ways EM field modes will be ignored in comparison to 
the cavity loss rate T. Thus we have r 3> 7. The case of 
non-negligible 7 is discussed by [H[ for the case of a two 
level atom. 

For density matrix elements in the basis |n; v) (where 
n = 0, 1, 2, . . . is the cavity photon number and v = 2, 1, 
specifies the upper, intermediate, and lowest atomic 
states) the master equation can be used to derive a set 
of coupled linear equations for the density matrix ele- 
ments Pnu-mfi in which the coefficients depend on the 
cavity frequency uj c , the detuning 5 between the cavity 
frequency and the average atomic transition frequency 
(5 = lu c — ujq), the half difference between lower and up- 
per transition frequencies (5, the coupling constants <?i, 
<72 and the cavity decay rate T. The coupled density ma- 
trix equations are set out in a suitable form in Appendix 
IaI The equations make up almost independent coupled 
sets of density matrix elements, if coupling only via the 
reversible processes associated with gi , g 2 is considered. 
In general we designate the (n, m) set for all n, m to in- 
clude the density matrix element p„o : m o (associated with 
the atom being in the lowest state |0)) and any density 
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matrix elements coupled to it via reversible processes as- 
sociated with the atom-cavity mode interaction terms. 
As we will see, relaxation due to cavity decay T couples 
a set (n, to) to other sets. 

The sets with (n, m ^ 2) consist of 9 elements — 



J n-22: 



ra-22;m-ll> Pn-22; m0> f ra-11; m-22' 
Tl;mO' Pn0;m^22) PnO;m^Tl an d PnO; mO 



-11; m-U> 

and may be represented as a 1 x 9 column matrix 



p(n,m) 




(35) 



This set is based on the three states of the form In — 2; 2), 
|n — 1; l) and |n; 0) (n 2) which are coupled via the re- 
versible processes associated with the atom-cavity mode 
interaction. However the terms arising from the irre- 
versible processes associated with the relaxation term 

Tap a COUple in p^i 2 ; m^T2> Pn=l2;m,l> ^n^T2:m+10' 
Pnl; m-12J Pnl; ml) Pnl;ro+l0) Pn+T0;m^T2' Pn+TO; ml an< ^ 

Pn+To- m+i o' wn i cn are the density matrix elements in the 
(n + 1, m + 1) set. Such irreversible processes are of the 
form \n + 1; v) — > \n;u) and these link in the three states 
|n— 1;2) |n; 1) and |n+l;0) on which are based the 
[n + 1, m + 1) set of density matrix elements. 

Special consideration is required for the density matrix 
elements associated with the state 1 0; 0) (where there 
are no cavity photons and the atom is in the lowest 
state). This state is not coupled to other states via 
the atom-cavity mode reversible interaction processes, 
and its only link with other states is via the irreversible 
1 1 : 0) — > |0;0) relaxation process. The density matrix 
equations associated with the |0; 0) state take on special 
forms. Similarly, special consideration is required for the 
density matrix elements associated with the states |0; 1) 
and |1;0) (where there are no cavity photons and the 
atom is in the intermediate state or there is one cavity 
photon and the atom is in the lowest state) . These states 
are coupled together but not coupled to other states via 
the atom-cavity mode reversible interaction processes. 
Their only link with other states is via the irreversible 
1 1 ; 1 ) — > 1 0; 1) and |2;0) — > |1;0) relaxation processes, 
and 1 1; 0> is linked to |0;0) via the |1;0) -> |0;0) relax- 
ation process. The density matrix equations associated 
with the |0; 1) and |1; 0) states also take on special forms. 
We still find that the density matrix elements break up 
into separate coupled sets associated with the reversible 
atom-cavity mode interaction processes. Following the 
general procedure for designating the sets via the photon 
quantum numbers n, to associated with the density 
matrix element p„o ; mO included in the set, the (0,0) set 
only contains the single density matrix element /?oo;00) 



the (0,1) set contains two elements poo; oi, Poo-, 10, with 
the (1,0) set also containing two elements poi; 00) PiO; oo- 
The (0,2) set contains three elements poo;02, Poo; n 
and poo;20j the (2,0) set also contains three elements 
P02; 00) Pii; oo an d P20;00- Next there is the (1, 1) set with 
four elements poi ; oi,Poi;io,PiO;Oi and pi 0; io- Finally, 
there are two sets each with six elements, the (1,2) 
set consisting of p ou 02 , Poi ; n, Poi; 20, Pio ; 02, Pio ; 11 
and jOio ; 20i an d the (2,1) set consisting of 
Po2 ; oi) P02; 10, Pii; 01, Pii; 10, P20; oiand £>20;io- In- 
stead of there being 9 coupled equations as for 
the (n,m) set (n, m ^ 2), there are 1,2,2,3,3,4,6 
and 6 coupled equations respectively for the 
(0,0), (0,1), (1,0), (0,2), (2,0), (1,1), (1,2) and (2,1) 
coupled sets. The complete sets of equations can be 
obtained from the master equation. 

In all cases the density matrix equations can be con- 
veniently be written in matrix form as 

d 

—p(n,m) = —iA(n,m)p(n,m)+iB(n,m)p(n+l,m+l). 
at 

Expressions for the column vectors p(n, to) and the ma- 
trices A(n,m) and B(n,m) which specify the coupling 
within the (n, to) set and the coupling to the (n+1, to+1) 
set respectively, are set out in Appendix lAl for (n, m 2) 
and for the special (0,0), (1,0), (1,1), (2,1) and (2,2) 
sets. 

Taking the Laplace transform of the equations for 
p(n, to) gives 

(s+iA(n, to)) p(n, m)—iB(n, to) p(n+l, to+1) = p(n, to, 0), 

where s is the Laplace variable. The Laplace transform 
A(s) of a function Ait) which is bounded for t ^ is 
defined as for complex s in the right half plane via 



A(s) 



dt exp(-si) A(t) 



Re s > . 



(38) 



The form of these equations indicates that if relaxation 
is also taken into account, certain collections of the 
p{n, to) form independent coupled sets, the set contain- 
ing a specific p(n, m) also includes pin ± 1, to ± 1), pin ± 
2, to ± 2), ... , pin ± k, to ± k), . . ., where k is a posi- 
tive integer such that n ± k 0,TO±fc ^ 0. Thus 
p(0, 0), p(l,l), p(2, 2), . . .are all coupled, as are p(l,0), 
p(2,l), p(3, 2), . . .etc. We see from the Laplace trans- 
form equations that if p(n > N, m > M, 0) are all zero, 
then for n > N and to > M (s + iA(n, m)) p(n,m) — 
iBin, to) pin + 1, to + 1) = 0, and hence a solution to the 
equations given by p(rt, to) = for n > N and to > M 
exists. Since linear first order equations have unique so- 
lutions, then pin, m, t) = for n > N and to > M. This 
feature is the consequence of the relaxation processes oc- 
curing in one direction, with only the \n + 1; v) — > |n; v) 
transition (and not its reverse) taking place. 

The solution to the density matrix equations at time t 
can be related to any initial density matrix in terms of 
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the matrices U n ,m:d,k{t) representing the evolution super- 
operator U(f) for the master equation as 

p(n, TO,*) Un,m::l,k(t) p(h k, 0), (39) 

l,k 

where the t7 nm:! j ) fc(*) satisfy the same matrix equations 
as p(n,m,t), but have an initial condition involving a 
unit matrix E 



(40) 



The dimensionality of the matrices f7n,m::i,Jfc(*) and 
E( n .m) depends on the sets (n,m) and (l,k). For 
example, if both (n, to) and (I, k) are such that 
(n, m 2) and (I, fc ^ 2) then the matrix is 9 x 9. 
On the other hand for (n,m 2) and (l,k) = (0,1) 
or (1,0) the matrix is 9 x 2, whilst if (l,k) = (0,0) 
the matrix is 9 x 1. Other cases follow similar lines. 
In full for (n,m ^ 2) and (l,k 2) the 81 matrix 

elements of U njTn: a,k(t) are U( n v),(m\)::(ip),(ka)(t), 

where (nv) = (n — 22, n — 11, nO), (mA) = 
(to - 22, rn^T l, mO), (Z/3) = (/ — 22, Z — 11, Z0) and 
(ka) = (fc — 22, k — 11, fcO). Other cases follow similar 
lines. 

For the evolution matrices we can also obtain equations 
for their Laplace transforms, and we have 

(s + iA(n, to)) U(n, to :: I, k) 

-iB(n,m)U(n + l,m + l ::l,k) = 5(„, TO ),(i,fc)S( n ,m)- 

(41) 

As we will see, Laplace transforms U(n,m :: I, k) occur in 
the final expression for the spectrum. Note that as the 
p(n, to) form independent coupled sets with relaxation 
only linking neighbouring sets (n, to) and (n ± 1, to ± 1), 
the only non-zero U n ^ m -.:i t k(t) are such that n — m = l — k. 

For the initial condition of interest the atom is in the 
upper state |2) and the single mode is in the vacuum 
state and the reservoir of continuum quasimodes in the 
vacuum state we have 



Pi = PsiPri 
psi = |0;2)(0;2| 

pm = n(i°x°i)*' 



(42) 
(43) 
(44) 



A 



and hence the only non-zero pin, to, 0) is for n = to = 2 



p(2,2,0) 



(45) 



For the discussion above regarding the irreversible link- 
ing of the coupled sets of density matrix elements, the 
unique solution for this initial condition is such that all 
p(n, to, t) are zero except p(2, 2, t), p(l, 1, t) and p{0, 0, t), 
which are obtained from their Laplace transforms 

p(2,2,s) = (s + a(2,2))-V(2,2,0) 
p(l,M) = ( S + 2^(l,l))- 1 ii3(l,l)p(2,2,,s) 
p(0,0,s) = (s + ^A(0,0))- 1 ^B(0,0)p(l,l,s). (46) 

This solution involves inverting a 9 x 9 and a 4 x 4 matrix. 
For this case of finite initial energy, the spectrum will 
be be that associated with a finite photon pulse in the 
cavity, but as we will see the spectrum will be different to 
that for spontaneous emission from an excited three-level 
cascade atom in free space. 

For the evolution matrices that are required for deter- 
mining the spectra we have the solutions 

[7(1,0 :: 1,0; s) = (s + iA{\, 0))" 1 
LT(n+l,7i::l,0;s) = (n > 1), (47) 



and 



C/(2,1::2,1; S ) = 
17(1,0:: 2,1; a) = 

U(n+l,n:: 2,1; s) = 



(s + tA(2,l))- 1 
(s + iAil^))- 1 
xiB(l,0) {s + iA(2,l)y 1 
(n > 1). (48) 



This solution corresponds to the previous result that an 
initial density matrix with only non-zero elements in the 
(TV, M) set or below cannot evolve into a density matrix 
with non-zero elements in sets such as (N + 1, M + 1), 
(N + 2, M + 2), . . . , due to the irreversible nature of the 
relaxation processes. Further details are in Appendix IBl 



B. Detector atom inside cavity — Case A 

To evaluate the spectrum for the case where the detec- 
tor atom is inside the high Q cavity (Case A) and with 
initial conditions of the atom excited and the cavity mode 
empty of photons we take 



V- = p*a 

= P>*^2Vn \n- 1; v) (n; v\ 



(49) 



V+ = fitf 

= p \fm \m; A) (m — 1; A| . (50) 



The two time correlation function is then 
Tr SR (piV+(h) Vt(ti)) 



mA nu 



(51) 
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where the transition operators S m \ ]nv (t) = |m;A) (n;v\ 
are Heisenberg operators at time t, and the trace is over 
both system and reservoir states. The two time correla- 
tion functions can be evaluated using the quantum re- 
gression theorem [l2|, |4(J 5l[, the result for which is 
given in terms of the evolution operator matrix elements 
Um\: nw.if); ka (t) associated with the density matrix equa- 
tions and considered as a function of the time differ- 
ence r = | tx — ^2 1 0, and density matrix elements 



PmX;w(ti,i) considered as a function of the smaller of 
the two times t\ and t 2 . The derivation of the expression 
for the spectrum for Case A is given in Appendix IB1 

The spectrum is given by 

S(u) = S 2 (u) + S 6 {u), (52) 

with 



S 2 {lo) = 2Re-p-VT\/l x {pio ; io(e')^lO;00::Oi;Oo(-Ki> + e) + j O l;10(e')C / 10;00::10:00(-i^ + e)} (53) 

\p\ 2 ~ ~ ~ 

S 6 (uj) = 2Re-p- [\/T/302;ll(e'){^lt / 10,00::02;0l(-^ + e) + \Zl#il,01::02; 01 + e) + V2f7 2 0,10::02;0l(-iw + ()} 

+ V2 (OQ2; 20(e'){^l C/ 10,00::02; lo(~^ + e) + VTf/ll,01::02; lo(~i^ + e) + V^C^O, 10: :02; lo(-£W + t)} 

+ V / T j Oii ; ll(e / ){VT?7l0,00::ll;0l(-iW + e) + \/l^n,oi::ii ; oi(-iw + e) + V^CaO.lQnll.Ol (- «*> + e)} 
W2p u . 2o(e'){"/LC/lo ) oo : :ii,io(-«w + e) + V^^ii,oi::ii,io(-^ + e) + V2U20 ,10::11; lo(~*^ + e)} 
+ V / T/O20; ll(e'){V / lL''i 0i 00::20;0l(-*^ + e) + VT^11,01::20; 01 (~ »W + ,10::20; 0l( — + e)} 

+V2 P20; 2o(e'){\ / l^io,oo::204o(-^ + e) + \/l^n,oi::20,io(-iw + e) + \/2^20 ,10::20,lo( — + c)}]- 

(54) 

I 



We see that in Case A, the spectrum involves Laplace 
transforms of the evolution operator from the E/i,o::i,o(t)j 
C^i,0::2,i(*) an d f7 2 ,i:.-2,i(i) matrices and from the p(2, 2, t) 
and p(l,l,t) density matrix elements. The spectrum is 
calculated from the equation (|52|) using the results from 
equations 1(15]). (|4T|) and (|4*5|) . The spectral variable w 
occurs in the U matrices associated with (1,0) and (2, 1) 
sets. From the equations determining these, it is easy 
to see that the spectral variable appears via the factor 
Auj = (lu — u> c ), so the spectra can be shown as functions 
of the detuning of the spectral variable from the cavity 
frequency. The calculation requires inverting 6x6 and 
2x2 matrices as a function of the spectral variable ui to 
give the U terms and inverting 9x9 and 4x4 matrices 
to give the p terms. 

C. Detector atom outside cavity — Case B (end 
emission) 

To evaluate the spectrum in the case where the detec- 
tor atom is outside the cavity to detect end emission and 
thus samples the continuum cavity quasimodes (Case B) 
and with the same initial conditions as before we now 
take 

V- = J dAp c (A)fx*(A)b(A), (55) 

and will assume that (J,* (A) is also a slowly varying func- 
tion of A. 



Starting from the formal solutions of Heisenberg equa- 
tions of motion for the Heisenberg operators b(A, t) 
and a(t) and making a Markoff approximation based on 
slowly varying factors pc(A), pf(A) andH^*(A), we can 
show that V- is the sum of a free field term V£ and a 
cavity term V_ , where 

V- = Vl + VS (56) 

V F = J dA Pc {A)p*(A) exp(-iAt)6(A) (57) 

V c = M*a, (58) 
and where 

M , = ^ A Pc(A)M*(A)^(A) 

J u c -A + ie K ' 

is an effective dipole coupling constant. Thus the cavity 
contribution V_ to the spectral quantity V- is the same 
as for Case A, apart from a constant of proportionality. 
The details are given in Appendix [Cl 

The expression for the two-time correlation function 
will involve free evolution continuum quasimode con- 
tributions of the form Tr SR V£ (t 2 ) V F (f x )Y 
cross terms involving the continuum and 
cavity quasimode contributions of the form 

Tr SR (piVf(t 2 ) VS{h)) ,Tr SR (fr£f(i 2 ) V F (h)) 
and purely cavity quasimode contributions of the 
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form Ttsr \J>iV+ {ti) V^(ti)^J. Since the continuum 
quasimodes are in the vacuum state it is not difficult 
to see using results such as 6(A)(|0) (0|)a = and 
(|0) (0|)a^(A) = 0, that the continuum quasimode term 
and the two cross terms all give zero, leaving behind 
only the cavity quasimode contribution 



Tr 



SR 



(piV+fa) V_(ti)) = Tr SR (piVfih) \> c (ti)) 



where 



V c = M* V™ \n - 1; v) {n; v\ 



(60) 



(61) 



v£ = mJ2V™ K a > ( m — i; A| . (62) 



m A 



The spectrum for Case B will therefore have the same 
form as that for Case A. 



D. Detector atom outside cavity — Case C (side 
emission) 

In the case where the detector atom is outside the 
cavity to detect side emission, and thus responds to the 
atomic transition operators (Case C), we have 

V- = R 2 d 2 + R\ dx (63) 

= r;J2 k i) (»; 2 I + R t E I"- °> ^ : l ( 64 ) 

n n 

V+ = R 2 d+ + R l d+ (65) 



where we have used the notation R\ 2 — » R 2 , i?oi ~~ > Ri- 
Since both R 2 and 172 are both proportional to the scalar 
product of the vector dipole matrix element (2\ d |1) 
between the upper and intermediate state with either 
the polarization unit vector for the spontaneous emission 
modes or the cavity quasimode, and similarly for R\ and 
,9i in regard to the vector dipole matrix element (1| d |0) 
between the intermediate and lower state, it follows that 



R2 
Ri 



92 

9i 



(67) 



Since we have chosen g\ and g 2 to be real we will make 
the same choice for Ri and R 2 from now on. 
Hence the two-time correlation function is 



Tr SR [piV+{t 2 ) V^{t 1 ) j 

R\ ^2 (S m2 - t ml(t 2 )S n i^ n2 (ti) 
m n 

m n 

+R 2 R± ^2 ^2 (Sm2;ml(t2)SnO;nl(tl) j 
m n 

+ i?.li?2^E (S m i- m o(t 2 )S n i- n2 (ti)\. (68) 



Note the four contributions to the spectrum, two involve 
just one atomic transition, the other two involve both 
transitions. The possibility of interference terms is clear. 
The derivation of the expression for the spectrum is given 
in Appendix [D] 

The spectrum is given by 



S(oj) = 2Re-^f[ f7o2;Oi::20;Oi(-iw + e)p 2 0;02(e / ) + U 02 - oi :: ii; i {-iu + e) Pll; 02 (e') 

+J7o2;01::02;Ol(— *^ + e ) P02;02(e')] 
+ 2 Re pf [ f7oi;00::01;Oo(-*^ + e)P01;Ol(e) + #01; 00::10; 00 + e) p w . l (e ) 

+{#01;00::02; lo(~^ + e ) + #11; 10::02: lo(— + e)} p 2; 11 (V) 

+{C/oi ; oo::ii; io(-iw + e) + Un- i 0: :ii; io(-iw + e)} pii; n(e') 

+{#01;00::20; lo(~^ + e ) + #11; 10::20: lo(~^ + e)} P20; 11 (c')] 
+2Re-^|# 1 [ {f/oi;00::02;Ol(-iw + e) + C/n ; 10::02; 01 + e)} P02: 02 (O 

+{#01;00::ll;0l(— + e) + #11; 10::11; 01 (— + e)} Pll; 02 (^0 
+{#01;00::20;0l(— + e ) + #11; 10::20:0l(-«^ + e)} 020; 02 (c')] 
+2 Re ^P 2 - [ C/o2;01::02;lo(-«W + e)p 2:ll(e') + #02; 01::11; lo(-^ + e) Pll; 11 (e') 

+#02; 01::20; lo(~»W + e) p 2 0; ll(e')]- (69) 
I 

We see that as in Case A (and Case B), the spectrum in- volves Laplace transforms of the evolution operator from 
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the J7i,o::2,i(i), ^2,i::2,i(i) and Uifi-.:i,o(t), and from the 
density matrix elements p(2, 2, i) and p(l, 1, i). However, 
in Case C different elements are involved, so although 
the spectra will be qualitatively similar in the two cases 
there will be quantitative differences. The spectrum is 
calculated from the equation (|69|) using the results from 
equations (146[) . (|47[) and (|4"5| . The spectral variable w 
occurs in the U matrices associated with (1, 0) and (2, 1) 
sets. From the equations determining these, it is easy 
to see that the spectral variable appears via the factor 
Alu = (u — w c ), so again the spectra can be shown as 
functions of the detuning of the spectral variable from 
the cavity frequency. 



IV. CASCADE ATOM SPECTRUM: 
NUMERICAL RESULTS 

Results for the spontaneous emission spectra from the 
three level cascade atom in a damped high-Q cavity are 
presented in figures I5HTT1 and figures [TMT41 In all cases 
the spectrum S(u>) is shown as a function of the detuning 
Ao; = uj — uj c of the spectral variable uj from the cavity 
frequency uj c . A logarithmic scale is used for S(uj) in 
order to highlight the numerous spectral peaks, and Auj 
is in units of g = g2- The spectrum for Case A(B) 
detector in cavity or end emission — is calculated from 
equation (f52"| and the spectrum for Case C — side emis- 
sion is calculated from equation ([55]) . In all cases, except 
figures [T31 and [T31 the values of the coupling constants 
are g\ = g% = 1-0 for numerical computation. For Case 
A(B) the detector atom-cavity mode coupling constant is 
/i = 1.0 (M = 1.0), whilst for Case C the detector atom- 
cascade atom coupling constants are R\ = i?2 = 1-0 (ex- 
cept in figures [T3l and [T4 | . The spectra are shown for var- 
ious choices of the cavity detuning 6, the half difference 
between the two atomic transition frequencies 6, and the 
cavity decay rate T. Strong coupling cases g ^ T and 
intermediate coupling cases g ~ T are displayed in the 
figures. 

Spectral calculations can also be based on a derivation 
of the spectra using density matrix elements defined in 
terms of dressed atom states, rather than the uncoupled 
states that we have used here. This can result in sim- 
pler expressions in the strong coupling regime, where line 
width factors are small compared to dressed atom level 
splittings. Analytical formulae for positions, heights and 
widths of spectral peaks can be often obtained in the 
strong coupling regime. However, for calculations cover- 
ing the weak, intermediate and strong coupling regimes 
using the dressed atom basis does not gain much advan- 
tage and for the present the uncoupled basis will be used. 

Figure [3] is for the case of resonance, with 6 = 6 = 0.0, 
where the two atomic transition frequencies are equal 
and the same as the cavity frequency. A strong coupling 
situation with T = 0.1 is presented for both Case A(B) 
and Case C. A symmetrical six peak spectrum is shown. 
For Case C the spectral intensity becomes very small 
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FIG. 3: Cascade atom in high-Q cavity. SE spectra S(u>) 
versus spectral detuning from cavity frequency Alu = (to — 
uj c ). The case of resonance is shown, where each transition is 
resonant with the cavity frequency, 8 = S = 0. The coupling 
constants are g\ = g% = 1 and the detector coupling constants 
are fi, = Ri = Rz = 1. The cavity decay is V = 0.1. The solid 
line is for Case A(B) — end emission — and the dashed line is 
for Case C — side emission. 



when the spectral frequency equals the cavity frequency, 
i.e. when Auj = 0. This effect is not seen for Case A(B) 
and indicates the presence of interference effects. 

Figures H] and [5] are also for the same resonance situa- 
tion as in figure[3j 6 = 6 = 0.0 for Case A(B) and Case C 
respectively. In these figures the cavity decay rate T takes 
on the values T = 1.0, 0.1 and 0.01, traversing the regime 
from intermediate coupling to very strong coupling. In 
the strong and very strong coupling cases T = 0.1 and 

0. 01 the six peaks are very clearly seen, whilst for inter- 
mediate coupling only two peaks remain due to the line 
broadening effect of the larger cavity decay. Again, for 
Case C (figure [5]) the spectral intensity tends to zero for 
Au = and T = 0.1 and 0.01. However, for T = 1.0, the 
spectral hole is suppressed and in fact the spectrum in 
Case C is very similar to that of Cases A(B) for r = 1.0 
as seen in figure |U We know that for a two-level system 
in the low-Q limit the cavity field operators can be adi- 
abatically eliminated in favour of the atomic operators, 

1. e. the cavity field 'follows' the atomic state. For this 
reason we might expect that spectra based on field oper- 
ators (case A) and spectra based on the atomic operators 
(case C) would become similar in this limit. 

Figures [5] and [7] for Cases A(B) and Case C respec- 
tively, apply to the situation where the two atomic tran- 
sition frequencies are the same 5 = 0.0, but where the 
cavity detuning 6 may be non-zero. Cases where 6 = 0.0 
(resonance), —1.0 and —2.0 are shown. A strong cou- 
pling situation with T = 0.1 applies in both figures. The 
effect of cavity detuning is that the six peak symmetrical 
spectra, applying for resonance, are replaced by asym- 
metrical spectra still with six peaks, but with the hint 
that other peaks may be hidden in the relatively broad 
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FIG. 4: Cascade atom in high-Q cavity. SE spectra S(u>) 
versus spectral detuning from cavity frequency Aui = (ui — 
uj c ). The case of resonance is shown, where each transition is 
resonant with the cavity frequency, 8 = 8 — 0. The coupling 
constants are gi = g% = 1 and the detector coupling constant 
is jLt = 1. The cavity decay V = 0.01 is the solid line, F = 0.1 is 
the dashed line and r = 1.0 is the chained line. The spectrum 
for Case A(B) — end emission — is shown. 
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FIG. 5: Cascade atom in high-Q cavity. SE spectra S(uj) 
versus spectral detuning from cavity frequency Aui = (ui — 
uj c ). The case of resonance is shown, where each transition is 
resonant with the cavity frequency, 8 = 8 — 0. The coupling 
constants are g\ = g% = 1 and the detector coupling constants 
are Tii = R2 — 1. The cavity decay V = 0.01 is the solid line, 
r = 0.1 is the dashed line and T = 1.0 is the chained line. 
The spectrum for Case C — side emission — is shown. 



features. Again the hole in the intensity for Case C near 
Auj = remains when the detuning is non-zero, though 
it is not as deep as in figureOfbr T = 0.1 (dashed curve). 

Figures [5] and M for Cases A(B) and Case C respec- 
tively, also apply to the situation where the two atomic 
transition frequencies are the same, 5 = 0.0, but where 
the cavity detuning S may be non-zero. Now a very 
strong coupling situation with T = 0.01 applies in both 
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FIG. 6: Cascade atom in high-Q cavity. SE spectra S(u>) 
versus spectral detuning from cavity frequency Auj = (uj—uj c ). 
The case of equal atomic transition frequencies is shown, but 
with the cavity frequency is detuned from the average of the 
two transition frequencies. 5 = but <5 may be non-zero. The 
coupling constants are g\ = 52 = 1 and the detector coupling 
constant is fJ, = 1. The cavity decay is T = 0.1. The cavity 
detuning 8 — is the solid line, 8 — — 1 is the dashed line, 
8 = —2 is the chained line. The spectrum for Case A(B) — end 
emission — is shown. 
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FIG. 7: Cascade atom in high-Q cavity. SE spectra S(uj) 
versus spectral detuning from cavity frequency Auj = (uj—oj c ). 
The case of equal atomic transition frequencies is shown, but 
with the cavity frequency is detuned from the average of the 
two transition frequencies. 8 = 0, but 8 may be non-zero. The 
coupling constants are gi = 52 = 1 and the detector coupling 
constants are Ri = R2 = 1. The cavity decay is F = 0.1. 
The cavity detuning 8 — is the solid line, 8 = — 1 is the 
dashed line, 8 = —2 is the chained line. The spectrum for 
Case C — side emission — is shown. 

figures, so that the spectral peaks will become narrower. 
Cases where 5 = 0.0 (resonance) and —1.0 are shown. 
It is now clearly shown that the effect of detuning is to 
replace the six peak symmetrical spectrum for the res- 
onance case, with eight peaks. We also note that in 
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FIG. 8: Cascade atom in high-Q cavity. SE spectra S(uj) 
versus spectral detuning from cavity frequency Auj = (a>— ui c )- 
The case of equal atomic transition frequencies is shown, but 
with the cavity frequency is detuned from the average of the 
two transition frequencies. 8 = 0, but 8 may be non-zero. The 
coupling constants are g\ = 52 = 1 and the detector coupling 
constant is fi = 1. The cavity decay is V — 0.01. The cavity 
detuning 8 = is the solid line, 8 = — 1 is the dashed line. 
The spectrum for Case A(B) — end emission — is shown. 
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FIG. 9: Cascade atom in high-Q cavity. SE spectra S(uj) 
versus spectral detuning from cavity frequency Aw = (ui—cj c ). 
The case of equal atomic transition frequencies is shown, but 
with the cavity frequency is detuned from the average of the 
two transition frequencies. 5 = 0, but 8 may be non-zero. 
The coupling constants are g\ — gi = 1 and the detector 
coupling constants are R\ = R2 = 1. The cavity decay is 
r = 0.01. The cavity detuning 8 = is the solid line, 8 = — 1 
is the dashed line. The spectrum for Case C — side emission — 
is shown. 



the detuned case one of the peaks displays a distinctly 
dispersion-like line shape. The spectral hole at Auj = 
for Case C now reaches an even lower intensity (10 -4 
units) than for figure [7] (10 -2 units), and is equally sharp 
for the detuned case as for resonance. 



FIG. 10: Cascade atom in high-Q cavity. SE spectra S(u>) 
versus spectral detuning from cavity frequency Auj = (uj—oj c ). 
The case of cavity detuning from the average of the two tran- 
sition frequencies 5 = — 1. The atomic transition frequencies 
may be different. The coupling constants are g\ = gi = 1 and 
the detector coupling constant is fj, — 1. The cavity decay 
is r = 0.01. The half atomic transition frequency difference 
8 — — 1 = +8 is the solid line, 8 = is the dashed line and 
8 = 1 = —8 is the chained line. For 8 = +8, —8 the cavity 
frequency coincides with the lower, upper atomic transition 
frequency respectively. The spectrum for Case A(B) — end 
emission — is shown. 
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FIG. 11: Cascade atom in high-Q cavity. SE spectra S(u>) 
versus spectral detuning from cavity frequency Auj = (uj — 
uj c ). The case of cavity detuning from the average of the 
two transition frequencies 8 = — 1. The atomic transition 
frequencies may be different. The coupling constants are gi = 
<?2 = 1 and the detector coupling constants are Ri = R% = 1. 
The cavity decay is V = 0.01. The half atomic transition 
frequency difference 8 = — 1 = +8 is the solid line, 8 = 
is the dashed line and 8 = 1 = — 8 is the chained line. For 
8 = +8, —8 the cavity frequency coincides with the lower, 
upper atomic transition frequency respectively. The spectrum 
for Case C — side emission — is shown. 



Figures ITD1 and [TT1 for Cases A(B) and Case C respec- 
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FIG. 12: The uncoupled atom-cavity mode energy levels and 
the dressed atom energy levels allowing for atom-cavity mode 
coupling. The situation shown is that of resonance and with 
equidistant cascade atom levels. 5 = 5 = 0. The uncoupled 
states are denoted |n; v), where n = 0, 1, 2, . . . gives the num- 
ber of photons in the cavity mode and v = 0, 1,2 lists the 
cascade atom states. The three lowest dressed atom multi- 
plets (singlet, doublet, triplet) are shown. Allowed transitions 
between the dressed atom states are indicated. 



tively, apply to the situation where the cavity frequency 
is detuned from the average of the atomic transition fre- 
quencies, and where these atomic transition frequencies 
may themselves be different. The cavity detuning is 
5 = —1.0 in both figures. A very strong coupling sit- 
uation with r = 0.01 applies in both figures so that 
the spectral peaks will become narrower. Cases where 
S = —1.0 = +6, S = 0.0 (both atomic transition fre- 
quencies equal) and S = +1.0 = —6 are shown. Note 
that for the case 5 = —5, the cavity frequency coincides 
with the upper atomic transition frequency (u) c = u>i), 
whilst the case <5 = +<5, the cavity frequency coincides 
with the lower atomic transition frequency (tu c — 102)- In 
both Case A(B) and Case C up to eight spectral peaks 
occur. For Case C (see figure ITTj) an interesting feature 
is that the spectral hole at Aui — is not present when 
the cavity frequency is resonant with the lower atomic 
transition frequency (S = +5), but is still present when 
the cavity frequency is resonant with the upper atomic 
transition frequency. 

The same general features for the present Case A(B) 
and Case C situation in terms of the numbers of spectral 
features were found by Zhou et al. [HI], but the presence 
of the spectral hole at Aw = for Case C (which they re- 
ferred to as the emission spectrum) and its relationship 
to quantum interference was not reported in their pa- 
per. Their lineshapes were associated only with the finite 
spectrometer bandwidth incorporated in their definition 
of the spectra and did not reflect cavity decay — which 
was assumed to be zero. 

The spectral features can be interpreted in terms of the 
dressed atom model Q ■ Figure [12] shows the uncoupled 
atom-cavity mode states and their splitting into dressed 



atom states for the three lowest dressed atom multiplets 
in the case of resonance 5 = S = 0. The lowest mul- 
tiplet is a singlet and is based on the |0; 0) uncoupled 
states. The next lowest multiplet is a doublet and is 
based on the |0; 1) and |1; 0) uncoupled states. The next 
lowest multiplet is a triplet and is based on the |0;2), 
|1; 1) and \2; 0) uncoupled states. These are of course the 
only states that could be occupied in the cascade atom 
decay process. The uncoupled levels are separated by an 
energy huj c = Hloq. For the resonance case the averages of 
dressed atom multiplet levels does not change from the 
uncoupled energies, and the triplet levels are equispaced. 
For the resonance case, the energies of the thr ee triplet 
states relative to the uncoupled energy are +hy/2g\ + g\, 
0, —h^/2g\ + g%, whilst those for the doublet state are 
+%i, — hgi and for the singlet state. The position 
for the spectral features is in accord with these dressed 
atom energies. For non-zero <5, 8 the dressed atom mul- 
tiplets do not have these symmetry features. The selec- 
tion rules for electric dipole processes only allow transi- 
tions between neighbouring dressed atom multiplets. It 
is therefore easy to see from figure [12] that there will be 
up to eight different possible transition frequencies, all 
clustered around lo = uj Ci and in non-zero <5, 6 situations 
all eight frequencies could be present. In the resonance 
case 5 = 6 = the symmetry feature results in the middle 
to bottom dressed atom transition frequencies coinciding 
with two of the upper to middle dressed atom transition 
frequencies, thereby reducing the number of distinct lines 
to six in this case. The fact that photons of the same fre- 
quency could have come from either an upper to middle 
or a middle to lower dressed atom transition, and that 
these processes cannot be separately measured, is likely 
to result in quantum interference effects. 

To demonstrate the role of both upper to middle and 
middle to lower dressed atom transitions in the inter- 
ference process, figures [13l and [14] show the effect for 
the resonance case of switching off the middle to lower 
dressed atom transitions. The figures show the spectrum 
for Case A(B) — end emission — as the solid lines and for 
Case C — side emission — as dashed lines. In these two fig- 
ures the two atomic transition frequencies are equal and 
the same as the cavity frequency S = S = 0.0, the cavity 
decay rate is T = 0.01, the value of the upper transition 
coupling constant is g-i = 1.0 = g. For Case A(B) the 
detector atom-cavity mode coupling constant is fi = 1.0 
(M = 1.0), whilst for Case C the detector atom-cascade 
atom upper transition coupling constant is R2 = 1.0. In 
figure [T3l the lower transition and detector atom coupling 
constants are g\ = R\ =0.1, for figure [T4l the values are 
gi = R\ = 0.01. The solid lines in previous figures [4] 
(Case A(B)) and [5] (Case C) show the spectra for the 
resonance case with gi = R\ — 1.0, where the lower 
transition and detector atom coupling constants are the 
same as those for the upper transition. The interference 
hole at Auj = for Case C is present in figure [5] is much 
less prominent in figure [T3] and has disappeared in fig- 
ure 1141 as the lower dressed atom transitions have been 
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FIG. 13: Cascade atom in high-Q cavity. SE spectra S(u>) 
versus spectral detuning from cavity frequency Aui = (ui — 
uj c ). The case of resonance is shown, where each transition is 
resonant with the cavity frequency, 6 = 8 — 0. The coupling 
constants are gi = 0.1, c?a = 1 and the detector coupling 
constants are Ri = 0.1, R2 = 1. The cavity decay is T = 0.01. 
The spectrum for Case A(B) — end emission — is shown as the 
solid line, that for Case C — side emission — is shown as the 
dashed line. 
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FIG. 14: Cascade atom in high-Q cavity. SE spectra S(u>) 
versus spectral detuning from cavity frequency Aui = (ui — 
ijj c ). The case of resonance is shown, where each transition is 
resonant with the cavity frequency, 5 = 8 — 0. The coupling 
constants are g\ = 0.01, g% = 1 and the detector coupling 
constants are _Ri = 0.01, R2 = 1. The cavity decay is V = 
0.01. The spectrum for Case A(B) — end emission — is shown 
as the solid line, that for Case C — side emission — is shown as 
the dashed line. 

essentially switched off. In fact there is now a spectral 
peak at Aui = for both Cases A(B) and C for the situ- 
ation in figure [T31 This shows that both upper to middle 
and middle to lower dressed atom transitions are needed 
to produce the destructive interference effect associated 
with the spectral hole. 

The situation when the cavity frequency is resonant 



with the lower atomic transition frequency and signif- 
icantly detuned from the upper atomic transition fre- 
quency is associated with a dressed atom level scheme 
in which the upper triplet is replaced by a doublet as- 
sociated with the splitting of the degenerate uncoupled 
states 1 1 ; 1 ) and |2;0) together with a separate dressed 
atom state which is essentially |0; 2). The doublet based 
on the 1 0; 1) and |1;0) uncoupled states is still present, 
as is the singlet based on |0;0). In this situation for 
Case C (side emission), the spectral hole at Aw = is 
absent. For the situation when the cavity frequency is 
resonant with the upper atomic transition frequency and 
significantly detuned from the lower atomic transition 
frequency is associated with a dressed atom level scheme 
in which the upper triplet is replaced by a doublet as- 
sociated with the splitting of the degenerate uncoupled 
states 1 1 ; 1 ) and |0;2) together with a separate dressed 
atom state which is essentially |2; 0). The doublet based 
on the |0; 1) and |1; 0) uncoupled states is now replaced by 
two singlets, each associated with these uncoupled states. 
The singlet based on |0; 0) is unchanged. In this situation 
for Case C (side emission), the spectral hole at Aui = 
is prominent. Why this is the case and what accounts for 
the general absence of the spectral hole at Aui — for 
the Case A(B) (end emission) and its usual occurrence in 
Case C (side emission) spectra is not yet clear and a more 
extensive analysis is needed to find the answer, possibly 
based on dressed atom states. For the strong coupling 
regime, the use of dressed atom states may enable tran- 
sition amplitudes for different processes — beginning with 
the same initial state with the atom in the upper state 
1 2), no photons present in any mode and the detector 
in state \A) and ending in specific final states with the 
detector in state \B), to be evaluated analytically using 
a perturbative approach for the weak detector atom and 
cavity decay transitions, but treating the strong coupling- 
processes involving the coupling constants g\ , gi without 
such an approximation. Having analytic expressions for 
the various quantum pathways would facilitate an un- 
derstanding of the interference effects that emerge in the 
spectra calculated numerically via the master equation- 
quantum regression theorem approach. 



V. CONCLUSION 

In this paper we have examined the spontaneous emis- 
sion spectrum from a cascade atom located inside a cav- 
ity. That spectrum was examined for three different 
physical situations regarding the location of an idealised 
detector atom. In case A the detector atom was located 
inside the cavity and directly influenced by the cavity 
field. In case B the detector atom was exposed to the 
field emerging from a partially transmitting cavity mir- 
ror. In case C the detector atom responded to sideways 
emission from the cascade atom inside the cavity. The 
cascade atom was initially in its uppermost level and the 
cavity mode empty of photons, and thus a system with 
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two basic excitations was studied. The spectrum was 
found from the weakly coupled detector atom response 
and involved in each case Laplace transforms of the den- 
sity matrix and matrix elements of the evolution oper- 
ator in a super-operator form. The spectral line-shapes 
only reflect cavity decay, since the weakly coupled de- 
tector atom spectrometer had zero bandwidth. Spectra 
have been presented for intermediate and strong coupling 
regime situations where both atomic transitions are res- 
onant with the cavity frequency, for cases of non-zero 
cavity detuning and for cases where the two atomic tran- 
sition frequencies differ. 

The spectra for Case A and Case B were found to be es- 
sentially the same. The spectral features for Cases B(A) 
and C were qualitatively similar, with six spectral peaks 
for resonance cases and eight for detuned cases. These 
general features of the spectra could be explained via the 
dressed atom model. However, Case B(A) and Case C 
spectra differed in detail. In particular, the spectra for 
Case C exhibited a deep spectral hole when the spectral 
frequency was equal to the cavity mode frequency, a fea- 
ture that persisted, for strong coupling, over a wide range 
of detuning and coupling conditions. However, this spec- 
tral hole was absent for Case A, a feature that cannot 
yet be explained. The spectral hole in Case C was par- 
ticularly prominent for resonance conditions, when both 
atomic transition frequencies were the same and reso- 
nant with the cavity frequency. In this resonance sit- 
uation, photons associated with the middle to bottom 
dressed atom transitions coincide in frequency with pho- 
tons associated with two of the upper to middle dressed 
atom transitions, and therefore the origin of a photon 
detected in the spectrometer cannot be identified. This 
suggests that an interference effect is occuring, and this 
is confirmed by showing that the hole disappears when 
only the upper to middle dressed atom transitions are 
present. However, why the interference hole is absent 



in Case B(A) spectra is still unclear. One clue is that 
the spectral hole in Case C spectra disappeared in fig- 
ure [11] when the atomic transition frequencies were quite 
different and the cavity frequency was resonant with the 
lower atomic transition frequency. This disappearance 
was sensitive to both the detunings having opposite sign 
and a magnitude close to g — g\ = 52. The explana- 
tion is expected to be found from an examination of the 
dressed states in the strong coupling limit. This analysis 
is currently being undertaken. 
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APPENDIX A: COUPLED DENSITY MATRIX 
SETS 

For the (n, m) set (n, m 2) we define the 1x9 column 
matrix p(n, m) of density matrix elements in the (n, m) 
set as 



p(n, m) 




(Al) 



and introduce suitable 9x9 matrices A(n,m) , B(n,m) to incorporate the couplings, detunings and decay terms 
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with p = n + m, and 



B{n,m) = —iT 
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For the (0, 0) set we have 

p(0,0) = [p 0;Oo] 

4(0,0) =w c (0) [1] 

B(0,0) =[0 -iT] 
For the (1,0) set we have 
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For the (2,1) set we have 
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For the (2, 2) set we have 
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A(2,2) = w c (0) 



0000000 
10000000 
01000000 
00100000 
00010000 
00001000 
00000100 
00000010 
1 



- 


-92 









92 





















— 92 — 


(s+sy 


-Iff 


-91 \/2 









92 



















-9l V2 


Of „'T 









() 


92 










() 


92 










+(<5 


+s)- 


fir 


-92 





9i v2 













92 









-92 




-ir 


-91 \/2 





91 


V2 













92 









-91 \/2 ( 


-5+5)-|ir 










91 V2 













91 \/2 








+2<5-ir 




92 























91 %/2 





-92 


-(-«+*)-■ §iT 


-91 V2 


. 




















91 \/2 





-91 \/2 


-i2T 






- -ir 






















" 













-ir-s/5 




































-irV3 





















5(2,2) 













— ' 


Ty/2 


































-ir 2 

































- 


rVe 

































-ir-s/5 




































-ir%/6 













. 






















-ir3. 







(A17) 



APPENDIX B: SPECTRUM FOR CASE A 

For the cascade case A situation we take 

= fx*J2V^\n-l;v)(n;v\ (Bl) 

TIV 

V+ = fj,a 1 

= fj,^2y/m\m;X){m-l;X\. (B2) 

The two time correlation function is then 
Tr SR (piV+{t 2 ) V-(h) 



mX nv 



(A18) 



r 



where the transition operators S m \ ]n „(t) = |m; A) (n;v\ 
are Heisenberg operators at time t, and the trace is over 
both system and reservoir states. The two time correla- 
tion functions can be evaluated using the quantum re- 
gression theorem [12, 0, Hl|, the results for which is 
given in terms of the evolution operator matrix elements 
U m X; nv.-.ip; ka (t) associated with the density matrix equa- 
tions and considered as a function of the time differ- 
ence r = \t± — t%\ J? 0, and density matrix elements 
Pm\-,w(ti,2) considered as a function of the smaller of 
the two times t\ and t%. 



(B3) 



We find using the result (|3Tj) for £ 2 = ti + T ^ t\ 

Tr SR (piV + (t 2 )V-(ti))=\v\ 2 Y,H 
and using (j52")l for t\ — t 2 + r ^ t 2 

Tr SR (piV+(t 2 ) V-(t{) 



lA nv 



mX nv 
ini nv 



V^v'n E^AjST^lA-feaiK^TvC 1 ") P*a;m<(*2), 



(B4) 



(B5) 
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where to obtain the last result we have made the interchange mX <-» nv. 
Substituting into the result for the spectrum we have 



S(u) = 



1 

ft 2 



JJ dhdhexpiuih - t 2 )Tr s (pi V+(t 2 ) V-fa) 



fi 2 

2 



S(u) = 



J I dt 1 dt 2 expiuj{t 1 -t2)Trs (piV+{t 2 ) 

I 12 /■OO POO 

S ( UJ ) = -TT^2^2^^ / ^ 2 / *l CX P* W (*l" t 2)^ mA ;^TA :: to;— At) Pkcc;nv{t2) 

n mX nv ka J J *' 

\ll\ 2 f°° f°° 

+ T2 _ X)X]^ / ™X) / dtl dt 2 expiu(h -t 2 )U^ =lx . mX ..^ =ll/ . k(X {T) Pnv , ka {t 1 ) 

n mX nv ka J ° " '*i 

I 12 roo poo 

-TTJ2J2^^^2 dt2 dT e M^T)U mX .— X ., ka — v { T )pka;nv(t 2 ) 

n mX nv ka Ja J 

|/i| 2 f°° Z" 00 

h mX nv >~ J J 



(B6) 



To obtain the last expression, the integration variables 
have been changed in the first region (t\ 5? t 2 ) to r, t 2 via 
the transformation t x = r + t 2 ,t 2 = t 2 (so the Jacobian 
equals +1) and in the second region (t 2 ^ ti) to t\,r via 
the transformation t\ — t\,t 2 = t + t\ (so the Jacobian 
equals +1). 

Since by taking the complex conjugate of the density 
matrix equations we have 

Pka;nv{t) = Pnv:,ka(t)* (B7) 
U m \-, m-l\::ka: n-lv( T ) ~~ ^m-lA; mX::n-lv; ka ( T ) ' 

(B8) 



it is not difficult to see that the second term in the last 
equation is just the complex conjugate of the first, thus 
proving that our result for the spectrum is real. 



For each contribution the double integrals factorise, 
each giving a Laplace transform — albeit for s on the 
imaginary axis (which may need to be written as a lim- 
iting process). We find that 



S ( W ) = ^YY^^Y U ra\;—X,k a ,—vi^ + e )Pka- nv{e') 
mX nv ka 

\P\ 2 ~ ~ 

+ n/W™ U ^1X; mXr^lv; ka(+ iuJ + e ) Pnv; ka (f?) 

mX nv ka 

\p\ 2 ~ 
= 2Re J^YY^^Y U ^^^ka:T^v(- i ^ + e )ha-,nv(e), (B9) 

mX nv ka 

I 



where the limits e, e' — > are understood. The last result 
follows from the second term being the complex conju- 
gate of the first. 

The initial conditions will lead to restrictions on the 
terms to be summed over. Of course the quantities a, A, v 
already only sum over the three atomic states. As shown 
previously for the present initial conditions, the only non 

zero Pka; nv 

(t) are those in the (2,2), (1,1) and (0,0) 



coupled sets — for (2,2) these arc /9 2 ; 02, Po2 ; n, Po2 ; 20, 

Pll;02, Pll; 11, Pll;20, P20;02, P20; 11 and p 20 - 2 0] for the 

(1,1) set we have Poi ; oi,Poi ; io,Pio ; oi and pi 0; io and for 
the (0, 0) set we have poo-, oo(t) — corresponding to the sys- 
tem states |0; 2) , |1; 1) and |0; 1) being the only ones pop- 
ulated during the decay process. This will restrict the 
sums over k,n, and the pka;nv{e') that are involved are 
all Laplace transforms of members of the (2, 2), (1, 1) and 
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(0,0) sets for the p(n,m,t). In fact given that we also 
require n ^ 1 due to the ^/n, factor, the only /Ofc a;m ,(e') 
that are possible are /5o2 ; n, P02 ; 20, Pii ; n, Pn ; 20, P2o ; n, 
?20:20, which are all in the (2,2) coupled set and poi;iO) 
Pio ; io which are in the (1,1) coupled set. This then 
restricts the U mX .j^ x .. ka .^ v {-iu) + e) to be of the 
form. U mX .— x .. 02 . 01 {-iuj + e), U mX .— A::02; 10 (-iw + 

£ J> U mX ,— A::ll;0lH W + e), ^ mA; — A::ll;10(-* W + C ). 
^mA;^TA : :20;0l(-^ + £ )' ^mA; ^TA :: 20; 10 H W + £ ) &Ild 
^mA;^TA: : 01;Oo(-^ + e ), U mX; ^TA::10; 00 H W + e ) be- 



cause of the fc, n restriction (note the n — lv in the sub- 
script). Since p 2;01, P02;10, Pll;01, Pll;10, P20;01j P20;10 

are all in the (2,1) set these are Laplace transforms 
of the J7„, m: :2,i(t)- The other terms poi ; oo, Pio ; oo are 
both in the (1,0) set these are Laplace transforms of the 
U n .m::i,o(t)- The spectrum can be broken up into two 
contributions, Se(co) which is associated with pk<x;nv{e') 
for the (2,2) coupled set and £2 (a/) which is associated 
with Pka;nv{e') for the (1, 1) coupled set. This gives 



S{u) 



2RC "^2- E E vW^E^A; — A::fea;— „(~»W + c) Pfca; ^(0 
m)l,An)l,i/ fea 

S 6 (w) + S 2 (a;) 



S 6 (w) = 2 R,C-^- XI ^^l t/ mA;^TA::02;0l(-^ + e )?02;ll(e') 



m^l,A 
|2 



+ 2RC -^2- E ^^2 t/ mA;^TA::02;10(-^ + e )P02;20(e') 



m>l,A 



+ 2RC ^2- E V^^mA^A.ll^lH^ + ^Pll;!!^') 



m>l.A 



+ 2RC ^2- E V^^mA— A::ll;10(-^ + e)pll;20(e') 



+ 2RC -^2- E V^VT[/ mA; — A::20;0lH W + e )P20;ll( e ') 



i>l,A 



+ 2RC T2- E V^^mA— A::20;10(-^ + e)P20;20(e') 
m>l,A 



S 2 (w) = 2Re^ X VWlE/, 



mi; m-lA::01; 00 



(-iu + e) p i;io(e') 



i>l,A 



+ 2RC T2- E V^^ /lC/ mA;^A::10;00(-^ + e )Pl0;10(e')- 



r 



(BIO) 



(Bll) 



The six Se(uj) terms involve U n ,m::2,i(— * w + e) However, for the first six Sq (u) terms we see from the 

and the two S 2 {w) terms involve U n , m: :i,o(-iu + Laplace transforms 

e). Then because the possible U n ,m::2,i(t) and 

U n . m ::i.o(t) must be such that n = m + 1, we 

see that the non-zero t/ mA; ^=T A::fcaiS =T 1/ (-»w + e) for 

the first six Se(uj) terms U n ^ m ;-2,\{—iu) + e) must 

all be Laplace transforms from the sets of matri- , . .,, ~, „ „ ,. . „. ~,„ , 

ces U lfi J(t), U 2 ,,,, A t), Us,,,At), r,, :: ,,;/K ... and (« + °» ^ :: 2 ' D " °) ^ 2 ' 1 : 

the the non-zero V mX ^ x ,, ka — lv {-i^ + e) for the (* + *A(2, 1)) 17(2, 1 : 

last two S 2 (cj) terms f n ,m::i,o(— + e ) mus t all (s + iA(3, 2)) 1/(3, 2 : 

be Laplace transforms from the sets of matrices [s + iA(4 3)) £7(4 3- 

^7l,0::l,o(*)) f72,l::l,o(*)) ^73,2::l,o(*)> ^4,3::!, (*) > ■ • ■ 



2,l)-iB(2,l)[/(3,2 : 
2, 1) - iB(3, 2) £/(4, 3 : 
2,1) -^(4,3)^7(5,4 : 



2,1) = 

2,1) = £ 6 
2, 1) = 
2,1) = 
(B12) 
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and it is easy to see that the solution to these equations 
is 

17(2,1:: 2,1; a) = (s + iA(2, l))" 1 

17(1,0 :: 2,1; a) = (s + iA(l, 0))" 1 

xiB(l,0) (s + ^(2,l))- 1 
U(n+ l,n :: 2,l;s) = (n > 1). (B13) 

This solution corresponds to the previous result that an 
initial density matrix with only non-zero elements in the 
(N, M) set or below cannot evolve into a density matrix 
with non-zero elements in sets such as (N + 1, M + 1), 
(N + 2,M + 2),... , due to the irreversible nature 
of the relaxation processes. This means that for the 
first six S 6 (u) terms the U mX .^[ X . :kct .^„(-iu + e ) 
must be from the f/2,i::2,i(t) and f/i,o::2,i(t) matrix 
elements only, and this places a restriction on m 
and A. Since the only density matrix element in 
the (1,0) set of the form p m \.^zr[\ is jOio ; oo then 
in fact the only terms are for m = 1 and A = 
(namely t/io,oo::02 ; oi(-"*> + e), f7io,oo::02 ; io(-«w + e), 

^10,00::ll;0l(-iW + e), ?7l0,00:: 11; 10 + e), 

^10,00::20;0l(-«W + e) and f/lO,00::20,lo(-«W + e)). 

Since the only density matrix elements in the (2, 1) set 
of the form p mX -^ix are pn ; oi and P20; 10 then the 
only terms are for m = 1 and A = 1 or m = 2 and A = 0, 
(namely U U ; 01-02; 01 (-"*> + e), ?7ii ; oi::02; io(-«w + e), 

^ll;01::ll;0l(-«W + e), ^7n ; oi::ll; w(-iu> + e), 

i/ii;0i::20;0i(— «w + e) and {7 11; £i :: 2o,io (— + e) for 
m = 1 and A = 1; then E/20,iO::02 ; oi(— + e), 

^20,10::02; w{~iu + e), ^20,10::!!; 01 + e), 



^20,10::ll;lo(-«^ + e), ^20,10::20; 01 + e) and 

^20,iO::20,io(^* w + e) for m = 2 and A = ). Thus each 
of the first six S%{w) terms produces three contributions. 

Also, for the last two 52 (w) terms we see from the 
Laplace transform equations 

(s + iA(l,0)) 17(1,0 : : 1,0) - iB(l,0) U{2, 1 :: 1,0) = E 2 
(s + iA(2,l)) C7(2,l : : 1,0) - iB(2, 1) 17(3,2 :: 1,0) = 
(s + L4(3,2))f7(3,2 : : 1, 0) - iB(3, 2) [7(4, 3 :: 1, 0) = 

(B14) 

and it is easy to see that the solution to these equations 
is 

(7(1,0 : : 1,0; s) = (s + iA(l, 0))" 1 
U(n + l,n : :l,0;s) = (n > 1), (B15) 

from the same considerations as for the first six Sq(oj) 
terms. This means that for the last two S 2 (oj) terms the 

U m \;-^l\:: ka — „H W + £ ) mUSt b ° fr ° m thc ^l,0::l,o(*) 

only, and this places a restriction on m and A. Since 
the only density matrix element in the (1,0) set of the 
form p m \.-^-ri\ is /Oio ; oo then in fact the only terms are 

for m = 1 and A = (namely £/io;00::Oi,oo( — ico + e) and 
?7iO;00::io,oo(^^ + e ))- Each of the two S 2 (ui) terms pro- 
duces one contribution only. 

We finally have for the spectrum 

S(u>) = S 2 (uj) + S 6 (lj) (B16) 



S 2 (uj) = 2Rc— VT\/T{/5io;io(e') c/ io;00: : oi;Oo(-^ + e) + Poi; io(e')^iO; 00::10;00(-«W + e)} 
S 6 (w) - 2Rc^ 

X [vT/002; ll(e'){ A /l^l0,00::02;0l( — «W + e) + VTC7n ,01::02;0l(— + e) 

,10::02;0l( — ^ + e)} 
+ V2/?02;20(e'){^1^10,00::02 ; lo(-«*> + e) + >/lf7ll,01::02;10(-«W + e) 

,10::02; lo( — *W + e )} 
Pll; ll(e'){^^10,00::ll;0l(-iw + e ) + V^ll,01::ll; 01 + e) 

,10::ll,0l( _ + e)} 

+\/2pii ; 2o(e'){V^^io,oo::ii,io(-iw + e) + v/Tf7ii,oi::ii,io(-iw + e) 
+\/2J7 2 o,io :: ii; io(-«w + e)} 

+\/lp20;ll(e'){^^10,00::20;0l(-iw + e) + \/l£7ll,01::20; 01 (-"*> + e) 
,10::20;0l(— + e)} 

+V2 

P20; 20 ,00::20,lo( — *W + e) + VT£7 U ,01::20,10(— + e) 

+\/2J7 2 o ,10::20,lo(-* w + e)}]. (B17) 
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APPENDIX C: SPECTRUM FOR CASE B 

For the cascade case B we have 

V- = J dA Pc (A)p*(A)b(A), (CI) 

where p*(A) is also a slowly varying function of A. 

The Heisenberg equations of motion for the Heisenberg 
operators b(A, t) and a(t) are 

jb{A,t) = -iAb(A,t)-iW*(A)a(t) (C2) 
^a(t) = -iuj c a{t) - ihgd_(t) 



■i j dApc(A) 



W(A)b(A,t), (C3) 



and we easily obtain the formal solution of the first equa- 
tion as 



b(A,t) = exp(-iAi) b(A, 0) 



-iW*(A) / dr exp(-iAT)a(t-r). 
Jo 



(C4) 



If there were no coupling terms, the free evolution for 
a(t) follows from the second equation to be a(t) = 
a(0) cxp(— iu c t). 

Hence we see that V_ is the sum of a free field term 
Vf and a cavity term VS , where 



v. = v F + vS 

V F = J dA pc (A) li* (A) exp(-*At)6(A) 

V c = dr ^-i J dA Pc {A)p*{A)W*(A) exp(-iAr) ) a(l - r). 



(C5) 
(C6) 

(C7) 



For the cavity contribution VEi the quantity 
F(t) = (-iJdAp c (A)p,*(A)W*(A)cxp(-iAT)) in- 
volves slowly varying factors pc(A), p*{A) andTy*(A), 
so it is not difficult to see that if the frequency width of 
the product of these factors is A^ then the overall quan- 
tity F(t) decreases to zero over a time scale tc ~ 1/Ac- 
We can then invoke the Markoff approximation on the 
basis that tq is small compared to evolution times for 
the system and approximate a(t — r) over the small time- 
scale tc over which the quantity F(t) is non-zero by its 
free evolution expression a(t — r) = a(t) exp(+iw c r). Ex- 
tending the r integral to oo and introducing the usual 
integrating factor exp(— er) gives the following Marko- 
vian result for VS 



V c = M*a, 



(C8) 



where 



M* = J 



dA 



pc(A)p*(A)W*(A) 
u> c — A + it 



(C9) 



is an effective dipolc coupling constant. Thus the cavity 
contribution V2 to the spectral quantity V_ is the same 
as for Case A, apart from a constant of proportionality. 



APPENDIX D: SPECTRUM FOR CASE C 

For the cascade Case C we have 
V- = +i?i?r (Dl) 

= r;J2 i n ; x > ( n > 2 i + R * E K °) ^ x \ ( D2 ) 

n n 

V+ = i? 2 £++i?i£+ (D3) 

m m 

(D4) 

where henceforth we will choose i?iand R 2 to be real. 
Hence the two-time correlation function is 



Tr SR (piV+(t 2 ) V_(ii) 

= Ri 

m n 

+ ^lEE (^S m l;mo(t2)S n O;nl(ti) 
m n 

m n 

+ i?li?2^^] (S m i- m0 (t2)S n i- n 2(ti)\. 



(D5) 



24 



For the first term using the quantum regression theorem 
result (|3ip we find that for t 2 = t± + r ^ t\ 



Ri 

m n 

; m2::nl; ka 



m n ka 



(D6) 



and using f32|) for = i 2 + t ^ i 2 

fij ^<S'm2;ml(^2)<5'nl;n2(^l) / 
m n 

= R\ J]] J]] J]] U n 2- nl-.-.ka; m\(j ) Pka; m2(^) 

m n fe Q 

= -R-2 /J /J /J C^m2; ml-.-.ka; nl( T ) Pfca; 712(^2)- 



m n ka 



(D7) 



For the second term using the quantum regression theo- 
rem we find that for £ 2 = t\ + r t\ 



^1 \Sml-,mo(~t2)SnO;nl(tl) 
m n 

; ml::nO; ka (r) Pnl; ka 



m n ka 



(D8) 



and for t\ = t 2 + t ^ i 2 

^1 (^ ml \ mo(t2)S n o- nl{tl) ^ 

m n 

; nO::A;a; mO (r) Pka; ml 

m 11 feet 

= fl iEEE^ 

; mO::fco(; nO Pka; nl 



m n fea 



(D9) 



For the third term using the quantum regression theorem 
we find that for i 2 = tj + t <i 



m n 

= R2R1 t/ml;m2::n0; fca( r ) Pnl; fea(*l)j 

(D10) 



m n fca 



and for ti = t% + r ^ i 2 



#2-^1 y»5'm2;ml(^2)<S'nO;nl(' ; l) 
m n 

= i? 2 i?l J7 n l ; „0::/ca; ml(f ) Pka; mife) 

m n ka 

= R2R1 ^ ^ ^ J7,„i (t) Pka;n2(t 2 )- 

m n ka 

(Dll) 

For the fourth term using the quantum regression theo- 
rem we find that for i 2 = t\ + r ^ <i 



RlR2 yS'ml; mo(i2)<5„i ; „ 2 (ii)^ 

m n 

= -R1-R2 ^ ^ ^2 ^ m °; 

ml::nl; &a Pnl; ka (tx), 

(D12) 



m n fcct 



and for ti = t 2 + r ^ i 2 



RlR% yS' m i :m o(i2)<5'nl;n2(il) / 
m n 

= R1R2 J7«2; nl::fcce; mo( T ) Pfcct; ml(*2) 

m n fca 

= R1R2 Pm2;ml::/:a; no( T ) P/ca; ill (£2)- 

(D13) 



m n fea 



For the spectrum we have 



S(w) 



/; 2 



1 

"ft 2 " 



dt^expicoit! -t 2 )Tr s [piV + (t 2 ) V-(h) 
dti«ft 2 expiw(ti - i 2 )Tr s (p* K + (t 2 ) V-(<i)) 



*2> tl 
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and substituting gives 

S(u) = S(u>) tl>t2 +S(u 

f£Z ^ ^ ^ roo roo 

S(u)) tl ^t 2 = TJ^UIX!/ dt2 / dtiexpiu(tl - £ 2 ) ^ m2 ;ml::to;r l l(T)/9to;n 2 (i 2 ) 
n m n ka J J ^ 



Thus 



j^2 ^ ^ ^ roo roo 

/ dt2 / dtl expiuj(ti - t 2 ) U m l; mOr.ka; nO (t) Pka-.nlih) 
n m n ka J ° Jt * 

^R 2 Ri 



R1R2 



poo roo 

- y^y^ y^ / ^ 2 / dtiexpiw(ti-t 2 ){7 m i ; mO::b; nl (t) pka;n2(h) 

m n ka ^° 
f£ jj> ^ ^ ^ roo roo 

'-^r~ y^.yi y^ / d< 2/ dt 1 G^iuj{t 1 -t 2 )u m2 - m i y ± a ., n Q(T)p ka , nl {t 2 ) 

71 m n ka J ° Jt * 
R 2 ■■ /■oo ^00 

%), 2>tl = d<1 / <&2expiw(ii -t2)f/ml;m2::nl;fca(T)Pn2;fca(*l) 

jj>2 ^ ^ ^ roo roo 

^2 ^2 ^2 / d<i / d<2 cxp ~~ t2 > Um0 

; ml::nO; fea \J~ ) Pnl; ka 

h m n ka J ° J *i 

m2::n0; fca (t) Pnl; ka (*l) 

ft m n to ^ 
j> jj> ^ ^ ^ roo roo 

■Tj^y^y^y] / ^1/ d*26xpiu;(ti-t2)y m 0;ml::nl;ka(^)Pn2;fca(*l)- 
h m n ka J ° J ^ 



R\R 

+ 



R 2 ■■ roo roo 

S (u)t 1 ^t 2 = / d<2 / dT CX P( iuJT ) U rn2- m l::ka;nl(T) Pka:n2(t 2 ) 

n m n ka J( > J ° 

j^2 ^ ^ ^ roo roo 

+ ^^2^2^2 / dt2 / dT ex P(* wr ) U m l;mO::ka:no(T) Pka;nl(t2 
h rn n ka J " J ° 



R2R1 
~ti 2 ~ 

R\R 2 



+ 



h 2 



^2^2 ^2 dt2 I dT ex -P( iuJT ) U ml ; mOr.ka; nl (t) Pka;n2(t 2 ) 

m n ka Jo J ° 

[•CO pOO 

^2^2 ^2 dt2 dT ex P( iuJT ) P m2; ml::fca; nO (r) Pfca;nl(*2) 

m n to - 70 



Jj>2 ^ ^ ^ rOC roo 

S(u) t2>H = Tfy^y^X!/ dtl / rfT ex PM wr ) U m l;m2::nl; to (r) p„2; to(*l) 
^ m n to 70 J ° 
j^2 ^ ^ ^ roo roo 

+ "p" ^2^2^2 dtl dT CX P(-^ r ) U raO ; ml::nO; fca (r) Pnl; ka 

(*i) 

ft m n to 7o - 70 

^ ^ ^ roc />oo 

^ — w~ ^2^2^2 dti dr cx pH wt ) u mi ; m2::n0; fca (r) Pnl; ka (ti) 

11 m n ka J J 

JllH 2 roo roo 

H r^— 5Zy^y^ / d< l / rfr CXp(-iwr) f7 m 0;ml::nl;fca(T")Pn2;fea(*l)- 

n m n ka J ° J 



(D14) 



(D15) 



(D16) 



To obtain the last expressions, the integration variables have been changed in the first region (ti > t 2 ) to r, t 2 via 
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the transformation t\ = r + t 2 ,t 2 = t 2 (so the Jacobian it is not difficult to see that the terms in the last equation 
equals +1) and in the second region (t 2 ^ ti) to t\,T via occur in complex conjugate pairs, thus proving that our 
the transformation t\ = t\,t 2 = r + t\ (so the Jacobian result for the spectrum is real, 
equals +1). 

Since by taking the complex conjugate of the density 
matrix equations we have 

Pka:nv{t) = Pnv;ka{t)* (D17) 
U m \ ; -^~r[\;;ka;7nv( T ) = U-^^ x . mX ..^l u . ka {T\T),l%) 



Hence we find that 

fOO />oo 



S(UJ) = 2 R e Tf EEE / dt2 / dT ex P( iuJT ) U m.2; m l-.:k a -,nl(T) Pk a -,n2(t2) 
11 m n ka J ° J ° 
j^2 ^ ^ ^ roo r-co 

+2Rc ~d"EEE / dt2 dr exp(iuJT)U ml:m0::ka . n0 (T) p ka . nl (t 2 ) 
n ™ „ t„ Jo Jo 



m n ka 
j£ j£ ^ ^ ^ roo i-ao 

+ 2Re ^EEE / dt2 / dT exp(iujT)U m i ;m O::ka;nl(r) Pka;n2(t2) 
n m n ka J ° J ° 



+ 2 RC 1 2 EEE / dt2 / dT exp(iujT)U m 2; m l::ka;no(T) Pka-.nlih)- 

n ™ „ Jo Jo 



m n ka 



(D19) 

For each contribution the double integrals factorise, each giving a Laplace transform albeit for s on the imaginary 



axis (which may need to be written as a limiting process). We find that 

fig 

h 2 

+ 2Re — j E E U m l;mO::ka;no{—iw + e) Pka;nl{c') 



R 2 ~ 

S(u) = 2Re TfEE '^2 U m2;ml::ka;nl{-i^ + e)pka;n2(e') 



m n ko 

?2 



h 2 

m n ka 



R R 

(-iuj + e) pka,n2{t ) 



h 2 

m n ka 



R R 

; ml::fca; nO 

(-i0J + e) pka;nl(e'), 



h 2 

m n ka 



(D20) 



where the limits e, e' — > are understood. 

The initial conditions will lead to restrictions on the 
terms to be summed over. Of course the quantity a 
already only sums over the three atomic states. As 
shown previously for the present initial conditions, the 
only pka:, n 2{t) and Pka;ni{t) that could be non-zero 
are those in the (2,2), (1,1) and (0,0) coupled sets — 
for (2,2) these are po2;02, P02 ; n, /Oo2 ; 20, Pn ; 02, Pii ; n, 

Pll;20, P20;02, P20; 11 and p 2 0:20', for the (1,1) Set WC 

have poi;Oi,Poi;io,Pio ; oi and pi 0; io and for the (0,0) set 
we have poo; oo (t) — corresponding to the system states 
|0; 2) , |1; 1) and |0; 1) being the only ones populated dur- 



ing the decay process. The only non-zero elements of the 
form pka;n2{t) arc /o 2:02, Pii:02 and p 2 0;02 from the (2, 2) 
set. The only non-zero elements of the form pka-.ni(t) 
are p 2 ; n, Pn ; n and p 20; n from the (2, 2) set and poi ; oi 
and pw-oi from the (1,1) set. No elements from the 
(0, 0) set are involved. Only the Laplace transforms of 
these elements will appear in the expression for the spec- 
trum, and this restricts the related values of k, a and n. 
The first and third terms in the last expression for the 
spectrum will thus give three different k, a and n con- 
tributions, whilst the second and fourth terms will give 
five. Thus we may write 



S(w) = S 22 {lo) + S\i(lo) + S 2 i{oj) + S 12 (oj) 



(D21) 
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R 2 ~ ~ 

#22 M = 2Rc "?^f X! [^m2;ml::20;0l(-iw + e) / O 2 0;02(e') + ^m2; ml::ll; 01 + e)pil;02(e') 

m 

+t/,n2;ml::02;0l(-iW + e) /O 02; 02(e')] (D22) 

R 2 ~ ~ 

Sn(uj) = 2Rc-^-^ [^ml;m0::O2;10(-«W + e) / O02;ll(e / ) + U ml; m 0::U; H){-iu + e) Pll; 11 (e') 
m 

+^ml;m0::20; lo( — *^ + e) P20: 11 (0 + ^n»l; m0::01; Oo( — *W + e) Poi;Ol(e') 
+i/ml;m0::10;00(-iw + e) p W - 1 (e')] (D23) 

R R ~ ~ 

S 21 (u,) = 2Re^£ ^ 

I;m0::02;0l(-«W + e ) P02;02(e') + ^ml; m0::ll; 01 (~ + e ) Pll;02(e') 

m 

+C/ml;m0::20;0l(-iw + e) P20;02(e')] (D24) 

^ _____ 

Sl 2 (w) = 2 Re ^ [fm2;ml::02;10(-iW + e) P02;ll(e') + U m 2; ml::ll; lo(-«W + e) Pll; 11 (e') 



+fm2;ml::20; lo(-«*> + e) /320; ll(e') 

+fm2;ml::01;00(-iw + e) p l;0l(e') + C^m2; ml::10; 00 (~ «W + e) PlO; 01 (e')] ■ 



(D25) 



This result can be rewritten to group separately the terms where the p are from the (2, 2) from those from the (1,1) 
set: 

S(u>) = S' 2 (uj) + S' 6 (lj) (D26) 
R 2 ~ ~ 

S' 2 ((jj) = 2Re-j ^ [^ r ml;m0::01;00(-iw + e)P01;0l(e') + U m i- m 0::W; Oo(-«W + e) p W - 1 (e')] 



ft 2 

+ 2 Re ^ ^ [t/ m 2;ml::01;00(-iw + <0 P01;0l(e') + ^m2; ml::10; Oo(-«W + e) PlO; 01 (e')] (D27) 



i? 2 ~ 

SgM = 2Re TjX! [ f/ m2;ml::20;0l(-«W + e)p20;02(e / ) 



ft 2 

m 

+fm2;ml::ll;0l( — * w + e ) Pll:0 2 (e') + t/m2; ml::02; 01 ( — * w + e ) P02;02(e )] 

R 2 ~ ~ 

+ 2 R ° ^2 \- Um1 ' ™0::02; lo(-"*> + e) P02; 11 (e') + f7 m l; m0::ll; lo(-»W + e) Pn ; n(e') 

m 

+f/ml;m0::20; 10 + e) ,020; 11 (e')] 

+ 2 Re ft2 ^ [f/ml;m0::02;0l(-«W + e) p 2;02(e') + U ml - m 0::U; 01 (~»W + e) /Oil; 02 (e') 

m 

+fml; m0::20; 01 + e) /?20; 02 (e')] 

+ 2 Rc — \- Um2 -' ™i:-m; io(-«w + e) p 2; n(e') 

m 

+fm2;ml::ll; lo(-«W + e) ftl; ll(e') + ^m2;ml::20; lo(-*W + e) p 20; ll(e')]. 



(D28) 



The ^(w) terms involve U of the form U ntm:: i t o since Poi;00 an d pio;00 are both in the (1,0) set. From 
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the discussion for the Case A spectrum the only non- 
zero f/„ ;m:: i ; o are of the form t/i,o::i,o- Since there 
are no elements in the (1,0) set of the form p m 2 ;m i 
then the second line terms f7 TO 2;mi::Oi;Oo( — + e ) an d 
U m 2;mi::W;Oo{—iu + e) arc both zero. As the only ele- 
ment in the (1,0) set of the form p m i ;m o is Poi:00 then 
the first line terms only involve one value of m, which is 
0. This gives 

R 2 ~ 

S' 2 (u) = 2Re-p-[Uoi;00::Oi ; oo(-iw + e)poi;Oi(e') 

+C/bl;00::10;00(-iw + e) PlO; Ol(e')]- 

(D29) 

The S' 6 (u) terms involve U of the form U n ,m::2,i since 

P20;01, Pll;01, P02;01, A)2: 10, Pll;10, P20; 10, P02;01, Pll;01, 

P20;0i, P02: io, Pii; io and ^20; io all in the (2, 1) set (three 
repeated terms). From the discussion for the Case A 
spectrum the only non-zero U n ,m-.:2,i are of the form 



£^2, 1-2,1 and J7i,o::2,i- Now in the (1,0) set there are 
no elements of the form p m 2-mi, and the only element 
of the form p m i:m0 is Poi;00- However in the (2,1) 
set there is one element of the form p m2;m i — which is 
Po2; oi, and there is also one element of the form 

Pml; mO 

which is io. Thus the first and fourth lines for S' 6 (iv) 
will only involve single U, namely ?7o2;Oi::20;Oi( — + e), 

^02;01::ll;0l(-|w + e) and U 2; 01::02j01 + e) in the 

first line, and f/o2 ; oi::02 ; io(-iw + e), J7o2 ; oi::ii;io(-iw + e) 
and f/o2;Oi::20; io(~ * w + e ) m the fourth line. How- 
ever, in the second and third lines two U are in- 
volved, namely i/oi;00::02;lo(-«*> + e), ^01;00::11; lo(-«W + 

e) and ?7oi;00::20;io(-iw + e) or f/n ; i 0:: 02 ; w{-iu + e), 
fii ; iO::ii;io(-«wJ- e) and Un- io::20; ioJ-iw + e) for the 
second line, and £/ i ; 00-02; oi (-iw+e) , ?7oi; oo::ii; oi (-«w+ 

e) and C/bl;00::20; 01 + e) Or J7ll;10::02;0l(-«W + e), 

f/ii;iO::ii ; oi(-iw + e) and f7ll;10::20;0l(-iW + e) for the 
third line. This gives 



2 



R 

S' 6 (U>) = 2Re-^| [t/()2;01::20;0l(-iW + e) ,020; 02 (e') + ^02; 01::11; 01 (~ «W + e) Pll;02(e') + ^02; 01::02; 01 (""*> + e) P02;02(e')] 

R 2 ~ ~ 

+ 2Rc-^- [{?7oi;00::02;lo(-iW + e) + ?7n ; 1 0:: 02; lo(-»W + e)} P02; 11 (e') 

+ {t/oi; 00::11; lo(-«*> + e) + Un-, i 0: :ll; lo(-»W + <0} Pll; ll(0 
+ {^01; 00::20; lo(~«W + e) + f/ll; 1 ::20; lo(-«*> + e)} p 2 0; ll(0] 

R R ~ ~ 

+ 2 Re ^ 1 [{C/bl;00::02;0l(-iW + e) + U U ; 10::02; 01 (~»W + e)} /9 2; 02 (e') 

+ {f/01;00::ll;0l(-«W + e) + Un ; 10::11; 01 (~»W + <0} Pll;02(e') 
+ {t/oi;00::20;Ol(-«W + e) + f7n ; 10::20; 01 ("«*> + e)} p 2 0;02(e')] 

R1R2 ~ ~ 
+ 2 RC ^ [f/02;01::02;10(-iW + e) ,002; 11 (e') + tfo; 01::11; lo(-»W + e) Pll;ll(e') 

+ ^02;01::20; lo(-»W + e) p 2 0; ll(e')]- ( D3 °) 

Combining the results for S' 2 (ui) and S' 6 (u) we obtain for the spectrum 

R 2 ~ ~ ~ 

S(U>) = 2 Re -j [i/02;01::20;0l(-iW + e) P20;02(e') + J/02;01::ll;0l(-iW + e) 011; 02 (e') + ^02; 01::02; 01 (-«*> + e) p 2; 02 (e')] 



ft 2 

R 2 ~ ~ 
+2Re [J7oi;00::Oi;Oo(-iw + e) poi ; oi(e') + J/oi;00::io ; oo(-iw + e) pio ; oi(e') 

+ {^01;00::02; lo(-«*> + e) + Un- 1 ::02; lo(~««> + e)} p 02; ll(e') 

+{foi ; oo::ii; io(-iw + e) + 10::11 . w {-iu + e)} pn ; n(e') 

+{^01;00::20; lo(-«*> + e) + f/ n . 10:: 20; lo(-«*> + e)} p 2 0; ll(e')] 

R2R1 ~ ~ 

+ 2 RC ^ [{t/oi;00::02;Ol(-«W + e) + Un ; 10::02; 01 (~»W + e)} p 2; 02(e') 

+{f01;00::ll;0l(-iW + e) + Ul\- 10::11; 0l(-«W + e)} j0ll;02(e') 

+{C^01;00::20;0l(-iw + e) + f/ll; 10::20; 0l(-«*> + e )} ?20; 02 (O] 
R^R 2 ~ ~ 
+ 2 R<3 fc 2 t^ 02 ; 01::02 ; 10M W + £ ) P02; 11 (e') + ^02; ■■]::'. I : lo(-^ + e) Pn ; H (e') 



ft 2 



+ t^02; 01::20: lo(-«W + e) p 2 0; ll(e')]- ( D31 ) 
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